Angular momentum nonconservation and conservation in quasiclassical Positronium 
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(Dated: December 13, 2010) 

It is shown that due to Thomas precession, angular momentum is not generally a constant of the 
motion in a quasiclassical model of the Positronium atom consisting of circular-orbiting point charges 
with intrinsic spin and associated magnetic moment. Despite absence of externally-applied torque, 
angular momentum is a constant of the motion only if the electron and positron intrinsic angular 
momentum vector components perpendicular to the orbital angular momentum are antiparallel and 
of equal magnitude. 

PACS numbers: 41.20.-q, 31.15Gy, 45.05.+X 



I. INTRODUCTION 

The timing and sequence of events in the development 
of quantum theory seem to be such that little attempt 
was made historically to incorporate electrodynamical ef- 
fects of the electron intrinsic spin into the classical atom 
model as pioneered by Rutherford [J 0] ■ In the absence 
of knowledge of the intrinsic spin and its magnitude in- 
volving the Planck constant (h), there was apparently no 
means to relate its occurrence as a proportionality con- 
stant in energy levels observed in atomic spectra, other 
than by postulate, as was done by Bohr Q to significant 
success. Subsequently, both Sommerfeld's |4j refinement 
and extension of the Bohr model, and de Broglie's [f| 
postulate of an electron wave character, also occurred 
prior to general recognition that the electron must pos- 
sess an intrinsic angular momentum with magnitude pro- 
portional to h. The published inference of the spin exis- 
tence by Goudsmit and Uhlenbeck Q in February, 1926, 
was nearly simultaneous with the inventions of Heisen- 
berg's matrix mechanics 0] and Schrodinger's wave me- 
chanics j8], in 1925 and 1926. Furthermore, although 
these first modern quantum theories were created with- 
out much if any knowledge of an intrinsic spin involving 
the Planck constant, it was soon incorporated in the new 
quantum theories by Pauli |9j, and its fundamental ne- 
cessity in quantum theory was recognized through the 
work of Dirac [10( published in 1928. 

Yet one must wonder what attempts might have been 
made to incorporate the intrinsic spin into the classical 
atom model of Rutherford, had not the modern quantum 
theory emerged so abruptly and successfully. The Bohr- 
Sommcrfeld model was itself remarkably successful, in its 
heyday, but had only a postulatory basis. Surely the dis- 
covery of intrinsic angular momentum and involving the 
Planck constant would have suggested the possibility of 
a dynamical basis for the Bohr postulate that orbital an- 
gular momentum is quantized in whole multiples of the 
reduced Planck constant, h = h/(2ir). It would have 
been natural to suppose that a basis for this postulate 
(as well as the two additional postulates involving h pro- 
posed by Sommerfeld soon after) might emerge from the 
classical electrodynamics of elementary particles carrying 



magnetic dipoles as well as electric charge. 

Instead, with the rapid emergence of modern quan- 
tum theory nearly simultanteously with the discovery 
that subatomic particles possess an intrinsic angular mo- 
mentum, interest in the old quantum theory based on 
modifications of the classical, point-charge Rutherford 
atom model rapidly waned. However, at least prior to 
Dirac's work, a strong impetus existed to look closely at 
the classical atom model with spin, due to the so-called 
spin-orbit coupling anomaly. The proposal of the spin by 
Goudsmit and Uhlenbeck had been based on atomic emis- 
sion spectrum measurements made by Zeeman for dis- 
charging gas under the influence of an externally-applied 
magnetic field. But, at the time of its introduction, the 
intrinsic spin idea had a serious problem. Electrodynam- 
ically, it did not account consistently and simultaneously 
for the magnitudes of the Zeeman effect and of the spin- 
orbit coupling that manifested as spectral line splitting 
even in the absence of externally-applied magnetic field. 
The spin magnitude inferred from the Zeeman effect is 
h/Air (i.e., h/2), while the calculation of the emission 
line splitting due to spin-orbit coupling, assuming this 
spin magnitude, obtained a result too large by a factor 
of two. 

It was L. H. Thomas who proposed [llj that the spin- 
orbit coupling anomaly could be accounted for by the 
relativistic effect now known as the Thomas precession. 
The Thomas precession reduces the rate of precession 
of the electron spin axis by a factor of about one-half, 
and assuming the classical-physics relationship between 
rate of precession and applied torque continues to apply, 
the spin-orbit coupling magnitude must also be approx- 
imately halved from its value expected in the absence 
of the Thomas precession. Thomas observed that this 
explanation required that "secular" , or orbit-averaged, 
total angular momentum be a constant of the motion, 
which he found to be the case. In recent decades how- 
ever it has been recognized that there must be included 
in the electrodynamics of magnetic dipoles an effect due 
to the presence of "hidden" momentum. When hidden 
momentum is included in the analysis, Thomas' result of 
secular angular momentum conservation can no longer 
be obtained for hydrogenic atoms with heavy nuclei that 
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lack a substantial magnetic moment [12J. 

Hidden momentum is mechanical momentum imparted 
to the charged current carriers of a current-loop magnetic 
dipole by an electric field [l3l - fl5j . It is equal and opposite 
the field momentum of the dipole magnetic field crossed 
with the electric field, so that the total of the field plus 
mechanical momentum of the stationary magnetic dipole 
in an electric field may vanish. Under any applied force, 
present an electric field, the equation of motion of a mag- 
netic dipole must account for the change in the hidden 
momentum that accompanies translaional and orienta- 
tional motion of the dipole, and so its omission leads to 
erroneous and nonsensical results. In particular, omit- 
ting the hidden momentum leads to violation of Newton's 
law of action and reaction in the interaction of magnetic 
dipoles with electric monopoles. Therefore, the quasicla- 
sical model of hydrogenic atoms considered by Thomas, 
where the electron intrinsic spin and magnetic moment 
is accounted for but the nuclear spin and intrinsic mag- 
netic moment is negligible, and the hidden momentum 
is omitted, is inconsistent. The attractive force on the 
electron toward the nucleus does not equal the attractive 
force on the nucleus toward the electron. 

Incorporating the hidden momentum does not entirely 
remedy the problems with the quasiclassical hydrogen 
model. Incorporating the hidden momentum in the elec- 
trodynamics obtains consistency between the two alter- 
native calculations of the binding force [H|-[l9[, but re- 
sults in nonconservation of the total angular momentum 
[2(i| . The rate of electron spin precession, which must in- 
clude the effect of the Thomas precession, in agreement 
with Thomas, is approximately halved compared to its 
value absent Thomas precession, but the orbit precession 
frequency is doubled compared to the value calculated 
by Thomas. The orbit precssion frequency magnitude is 
doubled due to the presence of the hidden momentum, 
and is unaffected by the Thomas precession that halves 
the spin precession frequency. Simple dynamical consid- 
erations then lead to an inescapable conclusion that the 
total angular momentum is not a constant of the mo- 
tion. The total angular momentum magnitude is con- 
stant but the vector total angular momentum precesses 
around a fixed direction, even absent an externally ap- 
plied magnetic field. This situation is not remedied ei- 
ther by accounting for the field momentum, nor by a 
careful relativistic treatment that takes proper account 
of the difference between the center of energy and center 
of mass. The field angular momentum can be shown to 
vanish identically for the hydrogenic atom [21( , and since 
the observation of angular momentum nonconservation 
can be made even for the minimally relativistic situation 
of quasiclassical Hydrogen with a circular-orbiting elec- 
tron at the ground-state Bohr radius, relativistic effects 
other than the Thomas precession are negligible and so 
cannot change this outcome. 

Beyond angular momentum nonconservation due to 
the Thomas precession, there is a problem of consistency 
of the spin-orbit interaction energy calculation, when the 



hidden momentum is included in the calculation, if it is 
assumed that the Thomas precssion causes a reduction of 
the spin-orbit coupling by the Thomas factor about one- 
half. The problem is that the orbital magnetic moment in 
the quasiclassical picture also precesses but isn't reduced 
by the Thomas precession, and so the energy associated 
with the orbit precession is not reduced by the Thomas 
factor, and so stands in disagreement with the calcula- 
tion focusing on the spin and including the Thomas fac- 
tor. This problem of inconsistency of description seems 
even more serious than the problem of angular momen- 
tum nonconservation and should lead to reconsideration 
of whether the Thomas precession can be accepted as 
plausibly reducing the spin-orbit coupling energy in the 
quasiclassical description. 

That the system consisting of a spinning magnetic elec- 
tron bound to a non-magnetic proton does not conserve 
angular momentum should perhaps not be considered 
surprising, in light of result obtained by Kiessling (22j, 
that classical electron theory does not conserve angular 
momentum in the absence of intrinsic spin. It therefore 
seems naive to expect it in the quasiclassical hydrogen 
atom model with a non-magnetic proton. Even though 
the proton possesses spin, without a magnetic moment it 
cannot couple into the electrodynamics as necessary to 
recover angular momentum conservation. The same con- 
sideration applies to accounting for the known intrinsic 
magnetic moment of the proton. As can be easily veri- 
fied, it is far too weak to correct the angular momentum 
nonconservation. 

The situation of angular momentum nonconservation 
in classical spin theory presents serious problems. One 
immediate consequence is that the picture that results is 
not even self-consistent, as already observed. There are 
two different spin-orbit coupling magnitudes that result, 
depending on the focus of the calculation. However there 
is another possibilty that might plausibly lead to a more 
satisfying result. In the hydrogenic atom model where 
the angular momentum nonconservation is observed, no 
accounting has been made of the compositcncss of the 
nucleus or proton. Considering even Hydrogen as the 
simplest case, and unknown to Thomas, it is now known 
that the proton consists of three particles with spin equal 
to the electron and presumably also possessing intrin- 
sic magnetic moment. It seems plausible that the dif- 
ficulties already described could be remedied if these 
could be taken into account, but the difficulties of do- 
ing such in the classical picture are considerable. This 
provides motivation to explore the quasiclassical electro- 
dynamical model of the Positronium atom. In Positro- 
nium, the problem of one particle lacking a useful in- 
trinsic magnetic moment is well remedied, yet the sim- 
plicity of the minimally-relativistic two-body problem is 
retained. However, one simplifying advantage of the Hy- 
drogen problem is lost in that there is no longer any ad- 
vantage in working in the electron rest frame. 

In the quasiclassical model of Hydrogen originated by 
Thomas and retained in many textbooks to the present 
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day, there is significant advantage to working in a ref- 
erence frame centered on the spinning particle, when 
the other particle can be considerd to be non-magnetic. 
Then, the more-complicated electrodynamics due to a 
moving magnetic dipole do not arise. This benefit is 
mooted in the analysis of Positronium, but nonetheless 
it is not too difficult to work directly in the laboratory 
frame. It can easily be verified that the techniques em- 
ployed including the motion of the magnetic dipoles yield 
the same result in a laboratory-frame analysis of Hydro- 
gen as is obtained in the electron rest frame. 

Based on the above motivating considerations, the ob- 
jective of the present analysis is to determine the angular 
momentum conservation properties, including any effects 
due to Thomas precession, of the quasiclassical Positro- 
nium model consisting of two point-charge particles, the 
electron and positron, also endowed with intrinsic spin 
and associated intrinsic magnetic dipole moment. This 
will be determined by developing an equation of motion 
of the total angular momentum that includes all terms 
to the order of the Thomas precession contribution. 

In the quasiclassical model, the particle positions and 
velocities, as well as their spin magnitudes and orienta- 
tions, are perfectly well defined at all times. To make the 
analysis as simple as possible, the particles are assumed 
to orbit each other circularly under Coulomb attraction. 
The orbital radius will be taken to be no smaller than the 
ground-state radius of the Bohr model of Positronium, 
which will ensure that relativistic effects other than the 
Thomas precession are negligible. Radiative effects and 
radiation reaction due to motions of the charges are also 
negligible over the course of a single orbit, as is the effect 
of delay. (These facts are demonstrated in an appendix.) 
It will then be readily seen, the total angular momentum 
is not generally a constant of the motion, for arbitrary 
orientations of the particle spins relative to each other 
and the orbital plane. Rather, angular momentum is a 
constant of the motion only if the particle intrinsic spin 
components in the orbital plane are aligned antiparallel 
and of equal magnitude. Thus, an expectation of angular 
momentum constancy gives rise to relative orientational 
constraints on the spins and orbit that may be akin to 
some behaviors of electron spins in atomic eigenstates in 
quantum theory. 

A sequel is planned to assess the system radiativity and 
radiation reaction forces related to the dipolar nature of 
the particles, and to determine if there exist forces that 
will make the special relative orientations preferred that 
obtain angular momentum constancy. It seems plausi- 
ble that in this way a new classical-physics mechanism of 
atomic radiative decay may be discovered, that may be 
compared with observation. However it does not seem 
likely that they will lead to general angular momentum 
conservation, because radiative fields are expected to be 
accompanied by radiation damping forces that them- 
selves balance the radiated angular momentum. Also, 
if these radiation fields due to intrinsic dipole motion are 
no larger than those due to the motion of the electric 



dipole due to charge separation, which are evaluated and 
shown to be inconsequential herein (see Appendices B 
and E), they can have no substantial effect. 

It should be stressed that the angular momentum con- 
servation violation according to the present analysis ex- 
ists despite accounting for angular momentum of the non- 
radiative electromagnetic field, which is not insignificant. 
Both the field and hidden angular momentum contribu- 
tions are essential in obtaining that the motion of the 
total angular momentum vanishes except for the contri- 
bution due to Thomas precession. 

Counter to the results presented herein, the expecta- 
tion that an analysis fully incorporating the Thomas pre- 
cession and all electromagnetic forces and fields should 
obtain conservation of angular momentum is bolstered 
by the work of Schild and Schlosser [23|, |24| , who showed 
that angular momentum is generally conserved in the 
electromagnetic two-body problem of charged particles 
with spin. Their manifestly-covariant general analysis ac- 
counts for both the Thomas precession and all (i.e. both 
nonradiative and radiative) electromagnetic fields. How- 
ever, the illustrative special case they considered, similar 
to the present system but restricted to aligned spins and 
orbit, was too restricted to uncover the interesting phe- 
nomonology to be described in the present work. This 
suggests a course of further investigation, to apply the 
theory of Schild and Schlosser to the situation investi- 
gated here. 



II. OVERVIEW OF THE ANALYSIS 

The analysis objective is to develop an equation of mo- 
tion for the total angular momentum of the quasiclassical 
Positronium atom model for the circular orbit and in- 
cluding the effect of the Thomas precession, for arbitrary 
relative orientations of the spins and orbit. The model 
is quasiclassical in the sense of being based on classical 
electrodynamics of point-charge particles with an ad hoc 
incorporation of intrinsic spin and magnetic moment. By 
assuming the quasiclassical Positronium atom is in the 
Bohr model ground state, the orbital velocity is small 
enough so that the primary significant relativistic effect 
is the Thomas precession. (There is also the relativistic 
effect that the moving magnetic dipoles acquire electric 
dipole moment as well.) It is shown in the appendix that 
any dynamical effects of radiation reaction and delay, as 
well as other relativistic effects, are negligible to the or- 
der of the present analysis. Therefore the central result 
of the present work, that according to the quasiclassi- 
cal analysis, total angular momentum is not generally 
conserved by electromagnetically-interacting point parti- 
cles with intrinsic spin, cannot be negated by a future 
analysis that includes radiation reaction and relativistic 
corrections to all orders. 

The total angular momentum to be considered has 
sources of four distinct types. These are the intrinsic 
spins, the kinetic orbital angular momentum due to the 
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motion of the particle masses, the so-called "hidden" or- 
bital angular momentum, and the field angular momen- 
tum. Equations of motion are developed for each source 
individually, and then summed to obtain an equation of 
motion of the total angular momentum. 

The equations of motions of the electron and positron 
intrinsic spins are developed in Section III using the 
equation of motion of the spin axis developed by 
Thomas [25]], which can be shown to be derivable 
from relativistically-covariant Bargmann-Michel-Telegdi 
(BMT) equation [26[. This equation of motion is 
straightforwardly specialized to the case of each particle 
moving in the electromagnetic field of the other. 

The equation of motion of the kinetic orbital angular 
momentum, that is, the orbital angular momentum of the 
particle masses, is developed in Section IV based on the 
various electrodynmic forces acting on the particles due 
to their attributes of electric charge and magnetic dipole 
moment. Importantly, due to the translational motion of 
the magnetic dipole carrying particles, they will acquire 
electric dipole moment that is not negligible to the analy- 
sis. Translational ( Stern- Gerlach) forces are then present 
due to the field gradients acting on the dipoles, as well 
as Coulomb and Biot-Savart forces acting on the charges. 
Additional non-negligible forces arise on the charges due 
to translational motions of the dipoles. The translational 
forces act to change the particles' mechanical momentum 
in accordance with Newton's law of motion, but the me- 
chanical momentum must include the hidden momentum. 
The time variation of the hidden momentum that must 
be accounted for in the equations of translational motion 
of the particles appears correspondingly in the equation 
of motion of the kinetic orbital angular momentum. 

When the general equation of motion of the particles' 
kinetic orbital angular momenta is available in terms of 
the total torque acting on the orbits, the next step is 
to identify all of the relevant translational forces and 
the torques they generate. After taking inventory of the 
forces and torques in Section V, they are evaluated ex- 
plicitly in Section VI. 

In Section VII the equation of motion of the total an- 
gular momentum is derived in two steps. The first step is 
to evaluate separately the rates of change of the kinetic, 
hidden, and field orbital angular momenta. These quan- 
tities are then summed to obtain an equation of motion 
of the total orbital angular momentum. The second and 
final step is then to sum the motion of the total orbital 
angular momentum with the motion of the electron and 
positron spin angular momenta to obtain the motion of 
the total angular momentum. The resulting equation of 
motion shows that the total angular momentum is not 
generally a constant of the motion, due to Thomas pre- 
cession, for arbitrary relative orientations of the particle 
spins and orbital angular momenta. Angular momentum 
is a constant of the motion only for certain specific rela- 
tive configurations of the spins and orbit. 

A set of appendices provide some standard results for 
reference and help to justify the approximations used. 



The analysis is performed in the laboratory frame 
throughout. 

A. Notation 

Particle positions relative to the center of mass are 
represented by r e and r p for the electron and proton 
respectively. In generic relations where the particle iden- 
tity is unimportant the subscript may be dropped for 
brevity, but in result equations r represents the vector 
from the positron to the electron, and r — \r\ to repre- 
sent the electron-positron separation. Thus r = r e — r p 
and for the equal-mass particles r e = —r p . Similarly, 
v = dr/dt will represent the electron velocity relative to 
the positron, while v e and v p will represent the particle 
velocities in the center-of-mass frame. 

Since there are orbital angular momenta from three 
contributing sources (due to motion of the masses, "hid- 
den" , and field) and due to each particle separately, it is 
important distinguish them clearly in the notation. The 
orbital angular momentum due to the motion of the elec- 
tron and positron masses, that is, the kinetic orbital an- 
gular momentum, is represented as L e and L p , and their 
sum L = L e + L p , is the total kinetic orbital angular 
momentum. The total hidden orbital angular momen- 
tum is represented as -Lhiddon = £hiddcn, e + -£<hiddon,p- 
The total field angular momentum is represented as 

-Afield = -Afield, e + -^field.p- 

B. Representation of Intrinsic Magnetic Moments 

The representation of the intrinsic spin and associated 
magnetic dipole assumed for the present analysis is an 
ideal classical current-loop magnetic dipole. The dipole- 
carrying particles are considered sufficiently separated so 
that only the magnetic far-field of the ideal dipole is sig- 
nificant. 

The intrinsic spin vector that follows from this repre- 
sentation has a well-defined orientation at all times. 

The magnitude of the spin angular momentum is as- 
sumed in the present analysis to be h/2. This is taken 
as the total magnitude, as opposed to being the magni- 
tude along a particular axis as in quantum theory. This 
magnitude assumption is primarily for convenience and 
the central result is insensitive to whether the spin mag- 
nitude is h/2 or (as in quantum theory) V3/2. 

It has also been suggested by various authors (see, for 
example [13] and references therein) that the electron 
intrinsic spin may be due to a classical but relativisitic 
circulatory motion of a charged point-like particle, and 
that evidence for this identification may be found in the 
so-called a ■ E term of the Dirac equation for a free elec- 
tron in an electric field. The existence of hidden mo- 
mentum and its assoctiation with intrinsic magnetic mo- 
ments of elementary particles, is not inconsistent with 
this classical "zitterbewegung" model for the spin. (This 
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use of the term "zitterbewegung" in this way should not 
be equated with the usage by Schrodinger in his analy- 
sis of the Dirac quantum mechanical relativistic electron 
theory.) The model used herein is explicitly assumed to 
have one property not possessed by the classical zitterbe- 
wegung model for the spin; that is, no time variation of 
the intrinsic fields that is faster than the classical orbital 
frequency of the particles. The effect of relaxing this as- 
sumption and considering a possibly more realistic model 
for the intrinsic fields (such as provided by Rivas (28[) is 
beyond the scope of the present work. However, as Ri- 
vas shows that the time-averaged fields of the relativistic 
classical zitterbewegung are asymptotic to the fields of 
an ideal current loop, it seems plausible that the results 
found herein may apply to the zitterbewegung model as 
well. 



C. Level of Approximation 

The level of approximation taken herein may be 
thought of as minimally relativistic in the sense that 
the relativistic 7 factor is usually approximated as unity. 
Two relativistic effects that are however central to the 
analysis and not neglected are the Thomas precession 
and the necessary electric dipole moment that accompa- 
nies a translating magnetic dipole. Attempts are made 
through analysis to justify neglecting other relativistic ef- 
fects such as propagation delay and the exact relativisitic 
description of the electromagnetic fields due to moving 
point particles, as the need arises. 

It is also assumed that the electron "g-factor," that is, 
the ratio of the electron magnetic dipole strength to that 
of a classical magnetic dipole constant current-carrier 
charge-to-mass ratio and equal angular momentum, is ex- 
actly two. That this assumption is only valid to a little 
better than the first order in v/c in the present applica- 
tion will be shown to be not significant to the result. 

The analysis is also restricted to circular orbit. This 
restriction is arbitrary but provides significant simplifi- 
cation. 



D. Explanatation of the Particular Sequence of the 
Analysis 

As stated above, the general analysis approach is to 
obtain separate equations of motion for each of the four 
identified angular momentum contributors, and then sum 
them to obtain the equation of motion of the total angu- 
lar momentum. It is to be expected however that these 
separate equations of motion of the angular momentum 
contributors are not independent of each other. For ex- 
ample and most obviously, the electromagnetic field an- 
gular momentum will depend entirely (and neglecting de- 
lay, which is defended in the appendix) on the instanta- 
neous kinematical description of the particles including 
their spin orientations and orientational motions as well 



as the particles' positions, velocities, and accelerations. 
The field configuration also affects the dynamical behav- 
ior through radiation damping forces, but fortunately ex- 
pressions for radiation damping forces are available from 
textbooks in terms of the charge and magnetic dipole ac- 
celerations, and in the current application these can be 
neglected, which is also shown in the appendix. Based on 
these considerations, the field angular momentum and its 
time rate of change will be evaluated as the last contribu- 
tion to the total, after the others are sufficiently defined 
so that the motion of the field angular momentum may 
be easily evaluated. 

The hidden angular momentum is determined by the 
kinematics and the instantaneous field configuration. 
The hidden angular momentum will be seen to have a 
description that is formally identical to the field angular 
momentum. Therefore it will also be evaluated after the 
kinematics have been specified by the equations of mo- 
tion of the spins and kinetic orbital angular momentum. 

Although the formula for the hidden angular momen- 
tum is formally similar to that for the field angular mo- 
mentum, hidden momentum differs from field momentum 
in an important respect. Unlike field momentum, hidden 
momentum enters directly into the equations of motion 
of intrinsically-magnetic particles, because any applied 
force must act on both kinetic and hidden momentum. 
In determining the equation of motion of the kinetic an- 
gular momentum, this will be handled by accounting for 
hidden momentum in the equations of translational mo- 
tion of the particles. 

An equation of motion of the spin orientations that de- 
pends only on the instantaneous fields is available from 
the literature. The orientation and motion of the spins 
will be seen to enter the dynamics directly and non- 
negligibly. Also, the hidden angular momentum and its 
motion will be seen to follow in a straightforward fash- 
ion from the spins orientations and motion, as well as the 
particle kinematics. Therefore the equations of motion of 
the spins will be evaluated as the first significant analysis 
step. Although the standard spin equation of motion is 
provided in terms of the fields, it will be convenient to 
rewrite them in terms of the other-particle relative po- 
sition, velocity, and spin axis orientation. The particle 
accelerations will not need to be included due to the neg- 
ligibility here of the acceleration field, as is shown in the 
appendix. 

Evaluation of the motion of the kinetic orbital angu- 
lar momentum follows sequentially after the motions of 
the electron and positron spins have been determined in 
terms of the kinematic description. When all the angu- 
lar momentum motion contributions are finally available, 
the motion of the total angular momentum is evaluated 
as their sum. 
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III. EQUATIONS OF MOTION OF THE SPINS derived by Jackson [29| from the BMT equation, and who 

also notes it to be equivalent to that developed earlier by 
The following equation of motion of the spin vector s, L. H. Thomas [HI]: 
for a particle with intrinsic spin and magnetic moment, is 



J 



ds 
~dt 



= — s x 



- 1 



B 



- 1 



-((3-B)(3- 



7 + 1 



/3 x E 



(1) 



where q is the particle charge, B is the magnetic field, 
E is the electric field, (3 — v/c where v is the particle 
velocity, and 7 = 1/ \/l — (3 2 , m is the particle mass, 
and t is time in the laboratory frame. This is Jackson's 
form except that the (signed) charge value of e has been 
replaced with q, so that e can be reserved for the funda- 
mental charge magnitude. Although this equation is de- 
veloped assuming an absence of field gradients, and there 
are clearly nonvanishing field gradients in the present ap- 
plication, the error of omitting the field gradient terms 
can be evaluated based on the work of Good [30j, and 
is not significant to the present analysis. In general, as 
will be shown explicitly, higher-order corrections to the 
motion of the spins are not significant to the present anal- 
ysis. It can be noted as well, Thomas applied Equation 
([1]) in his analysis reconciling spin-orbit coupling with the 
anomalous Zeeman effect, thereby disregarding the field 
gradients in the atomic domain. 

The equation of motion of the spin angular momentum 
vector of Equation (pj will be applied to determine the 
motion of the electron and positron spins in the quasi- 
classical model of positronium. Pursuant to this objec- 
tive it will be useful to first reduce the complexity of this 
equation in accordance with the approximations that are 
valid in the minimally-relativistic domain with respect 
to the quasiclassical Positronium atom in its Bohr model 
ground state. For example (see Appendix A), approx- 
imating the relativistic 7 factor as unity introduces an 
error of one part in 1//3 2 , with /? here approximately a 
factor of I /a 2 where a is the fine-structure constant, or 
about 5 x 10~ 5 . However, this level accuracy (that is, 
one part in a 2 ) is more than is needed here. Since the 
stated objective of the analysis is to obtain an equation of 
motion of the total angular momentum that includes all 
contributions significant to the order of the contribution 
of the Thomas precession, it's worth examining what is 
the contribution of the Thomas precession to Eq. ([1]) . 

The Thomas precession contribution to the motion of 
the spin vector in Eq. ((T|) is [cite Jackson?] 



-s x 



7 + 1 



(3x E 



-s x 



13 x E 



It will be seen below that in the present application 
(3 x E is of the same order of magnitude as B. Returning 



to consideration of Eq. ([T]), it can then be observed that 
the magnitude of the center term of the bracketed terms 
on the right hand side has a factor of /3 2 relative to the 
other two terms, and also relative to the contribution of 
the Thomas precession. Therefore this term is negligible 
in the present analysis, apart from that it will vanish 
under an approximation of the electron and positron g- 
factors as exactly 2. 

Approximating g as 2 in Eq. ([T]) is appealing simply 
on account of the algebraic simplifications it will provide, 
but the approximation that g — 2 is only accurate to 
slightly better (about one part in 10 3 here) than the first 
order in f3. An examination of the potential effects of tak- 
ing the approximation g = 2 is therefore warranted. Now, 
in the equation of motion of the total angular momen- 
tum, it is to be expected that the motion of the spins will 
enter directly in accordance with Eq. ([1} , since the total 
angular momentum itself will include the electron and 
positron spins directly. In this contribution the approxi- 
mation that g — 2 will have no untoward effect beyond a 
slight error of the net motion of the spins, that is negli- 
gible compared to the effect of the Thomas precession in 
Eq. ([T]). Additionally, the rate of change of orientation 
of the spin may enter also into the part of the motion of 
the total angular momentum that is due to the motion of 
the orbital angular momentum. However it will be seen 
in Section VII that the motion of the spins enters into 
the equation of motion of the orbital angular momentum 
only at order /3 2 compared to the instantaneous orienta- 
tions of the spins. Therefore, in the equation of motion 
of the total angular momentum that contains all terms 
to the order of the Thomas precession contribution, the 
contribution of the rate of change of the spins in the mo- 
tion of the orbital angular momentum is negligible. (This 
is shown explicitly in Section VIIB, Eq. (jl08l) .) There- 
fore there will be no further error introduced by taking 
the further approximation of the electron and positron 
g-factors to be exactly two, as opposed to their exact 
values of about one part in 10 3 larger, other than chang- 
ing the rate of motion of the spins slightly. This error will 
have no effect on the central result of the present analy- 
sis, that the total angular momentum is not generally a 
constant of the motion, as a consequence of the Thomas 
precession. Later it will be straightforward if desired to 
repeat the analysis with more precise g-factor values, to 
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confirm there is no effect on the central conclusion, but 
for the initial analysis it is a significant simplification to 
let 5 = 2. 

As justified by the considerations above, neglecting 
terms of relative order of /3 2 , and approximating 7 as 
unity and the g-factor as two results in the following 
equation of motion for the electron spin s e in the qua- 
siclassical Positronium atom, where E and B are the 
electric field and magnetic field of the positron: 



ds e 
dt 



mc 



B 



(2) 



where the generic q for the charge has been replaced 
with an electron-specific quantity q e so that the polari- 
ties of the electron and positron may be more easily ac- 
counted for. That is, q e = —e where e is the fundamental 
charge magnitude, and q p = e, where q p is the positron 
charge. Here and henceforth, m = m e = m p is the elec- 
tron or positron mass. 

The magnetic field at the electron is due to the in- 
trinsic magnetic moment of the positron and to the mo- 
tion of the positron charge. This can be expressed as 
B = Bp + B v . With the intrinsic magnetic moments 
modeled as classical current loop dipoles as described 
above, the magnetic field component due to the positron 
intrinsic magnetic moment is 



B2, this expression represents the exact magnetic field of 
the moving positron charge accurately to the first order 
in f3. Therefore by the same reasoning as above, this will 
introduce only an order /3 3 error in the motion of the 
total angular momentum.) Noting L = L 
obtains that 



B v = -%^L, (5) 

mcr a 

The electric field at the electron consists of the fieled 
due to the positron's electric charge plus an additional 
contribution due to the motion of the positron intrin- 
sic magnetic moment. The electric field due to transla- 
tional and orientational motions of the intrinsic magnetic 
dipoles is considered in detail below. The latter total 
magnitude is of order /3 2 compared to the Coulomb field. 
In general, the term in Eq. (|TJ) involving the electric 
field is of the same order as the term involving the mag- 
netic field directly, and therefore by the considerations 
stated above, that relative order f3 2 terms in the motion 
of the spin are irrelevant to the present analysis, it must 
be concluded that the contribution to the motion of the 
spin axes of the electric field due to motion of the intrin- 
sic magnetic dipoles is negligible to the present analysis 
of the motion of the total angular momentum. 

The Coulomb field of the positron at the electron is 



+ Lp — 2L p 



B M = 



3n (n ■ n 



(3) 



where fi p = q p s p /mc is the positron intrinsic magnetic 
moment, and n — r/r is a, unit vector in the direction 
from the positron to the electron. 

The magnetic field at the electron due to the orbital 
motion of the positron charge may be approximated as 



Qp 



2q p 



B v = -^v p x r = ^v p xr p = 



2q p 



3 u p 



p 3 V ' 

mcr a 



(4) 



with r = r e — r p = —2r p = 2r e , and L p — r p x mv p is 
the positron kinetic orbital angular momentum. (This 
expression for the magnetic field due to a translating 
point charge may be obtained from the exact Lienard- 
Wiechert magnetic field of a moving point charge by ne- 
glecting the acceleration (radiation) field, neglecting de- 
lay, and approximating 7 as unity. As shown in Appendix 



E 



(6) 



It is shown in the appendix that the description of the 
electric field as the instantaneous Coulomb field devia- 
tiates only by an order /3 4 error from the exact delayed 
Lienard-Wiechert velocity electric field of the electron or 
positron in the Bohr model ground state. It is also shown, 
the Lienard-Wiechert acceleration electric field term is of 
order /3 3 here and so is also negligible. Therefore the elec- 
tric field will be represented in the equations of motion 
of the spins as simply the undelayed Coulomb field as 
represented by Eq. ([6J , and the equivalent expression for 
the electron Coulomb field in the equation of motion of 
the positron spin. 

Based on the above considerations the equation of mo- 
tion of the electron spin vector of Eq. <j2j) becomes 



ds e 
~dt 



'3n (n- fi p ) -H p 



(7) 



where (3 e = v e /c is the electron velocity in the laboratory frame, and r = 2r e is the vector from the positron to 
the electron. 
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The total kinetic orbital angular momentum L is twice the electron kinetic orbital angular momentum L e 
r e x mc(3 e = —mc/3 e x r e , and using again that r = 2r e , so that mc/3 e x r = 2mcf3 e x r e = —2L e = L obtains 



ds e 
dt 



mc 



q P ( 3n (n • s p ) - s p 
mc V r 3 



q P 

mcr 3 



1p 
mcr 3 



(8) 



Replacing q p with its value of the fundamental charge 
+e, and q e by — e, and collecting similar factors obtains 
finally 



:S K x 



3n (n • s p ) — s ? 



(9) 



It is obtained similarly for the positron that 



p m 2 c 2 r 3 p 



Sr> X 



3n (n • s e ) — s e 



I" 



(10) 



It is also of interest to make note of what is the contri- 
bution of the Thomas precession to Equations (j9]) and 
(flTJ|) . The contribution of the Thomas precession to 
Equation is 



«TP,e 



r S R X 



5" 



(11) 



The contribution of the Thomas precession to Equa- 
tion (TO]) is 



brackets (and as modified by the Thorns precession when 
viewed in the laboratory frame). Since radiation inten- 
sity is proportional to the square of the magnitude of 
dipole acceleration, and for harmonic motion the accel- 
eration is proportional to the square of the frequency, 
this represents a very great increase in radiativity com- 
pared to the radiativity associated with the Larmor pre- 
cession of the spin-orbit coupling. Also, as is well known 
and will be treated explicitly in Section VI and else- 
where herein, the magnetic-dipole-carrying particles will 
acquire electric dipole moments due to their motion, and 
so acceleration of the spin vectors will potentially cause 
electric dipole radiation as well as magnetic dipole ra- 
diation. Such radiative mechanisms don't seem to have 
been considered previously and warrant further exam- 
ination to determine the degree of correspondence, if 
any, with observation. Dipole radiation due to motion 
of point charges is after all a good match for observation 
in the limit of large quantum numbers, as expected in 
accordance with the Bohr correspondence principle. The 
analysis here hints that there may be additional classical 
phenomena that can extend the correspondence further 
into the quantum domain. 



2 2„3 P 

mcr' 



(12) 



In Section VIII, Equations (J9J) and (flOf are used to 
determine the contributions of the electron and positron 
spin motions to the total angular momentum. An upper 
limit on the magnitude of the time rate of change of the 
spin vectors in the Bohr model ground state of Positro- 
nium is determined in Appendix D. This limit value is 
used to show the neglibility of certain effects and so help 
to simplify the analysis. 

Before proceeding, and although it won't be pursued 
further herein, it is perhaps worthwhile to make an ob- 
servation about the nature of the motion of the spins in 
the quasiclassical Positronium atom model, as embodied 
in Eqs. © and (TTU)) . It is interesting that there is an 
additional component in the motion of the spin vectors 
that is a much more rapid motion than the regular pre- 
cessional motions, that are essentially Larmor precession 
when viewed in the particle rest frame. Specifically, the 
first term in the square brackets of either equation varies 
at the orbital frequency, which, in the present applica- 
tion, is about eight decimal orders of magnitude larger 
than the Larmor precession frequency of the spin-orbit 
coupling, represented by the third term in the square 



IV. EQUATION OF MOTION OF THE KINETIC 
ORBITAL ANGULAR MOMENTUM 



In this section the equations governing the time rates 
of change of the electron and positron kinetic orbital an- 
gular momentum vectors are developed. As previously 
outlined in Section IID, the magnitude of the spin angu- 
lar momenta and the contribution of the Thomas preces- 
sion for the individual particles determines the accuracy 
needed for the hidden and field angular momentum con- 
tributions to the total angular momentum. In determin- 
ing the effect of the Thomas precession, there is no need 
to be concerned with terms that are orders of magnitude 
in /3 smaller than its own contribution. Also, although 
the kinetic orbital angular momentum is distinct from 
the hidden angular momentum, the presence of hidden 
momentum enters directly into the equation of motion 
of the kinetic orbital angular momentum, in a significant 
fashion. 

In the quasiclassical Positronium model viewed in the 
laboratory frame, the electron and positron are orbit- 
ing each other influenced by electric and magnetic fields 
due to the other-particle electric charge and dipole mo- 
ment. The predominant influence is the Coulomb attrac- 
tive force between the particles, but non-Coulomb forces 
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are present that cause deviation from Keplerian motion. 
Herein as noted it is assumed that the orbits are circular, 
but it cannot be expected that an exact circular orbit will 
result dynamically if the forces are not purely radial and 
conservative. However, assuming an initial condition of 
perfect circularity, the deviations from circularity on the 
time scale of an orbit can be shown to be small enough 
so as to be insignificant to the present analysis. 

The approach to determining the motion of the kinetic 
orbital angular momentum is as follows. First, the ba- 
sic equation of motion of the kinetic orbital angular mo- 
mentum is developed in terms of the kinetic linear mo- 
mentum that in the nonrelativistic regime is simply the 
mass times the particle velocity, separately for each parti- 
cle. This reduces the problem of determining the motion 
of the kinetic orbital angular momentum to one of de- 
termining the change of particle velocity in response to 
applied forces. Then, the equations of translational mo- 
tion of the particles are developed based on the applied 
forces acting in accordance with Newton's law of motion, 
where the momentum changing under the applied force 
is the total mechanical momentum that includes hidden 
momentum. Finally, the equation of motion of the total 
kinetic orbital angular momentum is developed in terms 
of the total torques acting on the individual particle or- 
bits. 



A. Motion of the Kinetic Orbital Angular 
Momentum of Each Particle Separately 

The positron kinetic orbital angular momentum is de- 
fined as 



L p = r p x P p — r p x mv p , (13) 

where P p is the positron kinetic momentum given in 
the nonrelativistic limit as P p — mv p . The motion of 
the positron kinetic orbital angular momentum is then 
given by 

L p = r v x mVp + r v x mbp — r p x mv p , (14) 
with v p = Tp. 

Similarly, for the electron, L e = r e x mv e . 

The equations for the motion of the kinetic orbital 
angular momentum of the individual particles show the 
need for evaluating the particles' acceleration. 



B. Hidden Linear Momentum and the Equation of 
Translational Motion of a Magnetic Dipole 

In general a force F on a classical rigid body is equal 
to its time rate of change of mechanical momentum, 
dPmcch/dt, where -P mcc h is the mechanical momentum 



of the body. In the nonrelativistic limit, and for a clas- 
sical body without a magnetic moment, the mechanical 
momentum is the kinetic momentum P = mv. How- 
ever, as described in the introduction, it is now generally 
recognized that if the body carries a net current in its 
rest frame, and if an electric field is present, then there 
is mechanical momentum in addition to the kinetic mo- 
mentum. This momentum is nonvanishing in the rest 
frame of the current-carrying body and is often termed 
"hidden" momentum. For a classical current loop mag- 
netic dipole of moment m, in an electric field E, the 
hidden momentum is Phidden = m x E/c [3l|, [32J. The 
existence of hidden momentum implies that the equa- 
tion of translational motion of a classical current-loop 
magnetic dipole must take into account the time rate of 
change of the hidden as well as the kinetic momentum of 
the body. Omission of the hidden momentum from the 
equation of motion will lead to a nonphysical dynamical 
description of the interaction of electric monopoles and 
magnetic dipoles. For example, if the hidden momentum 
is not taken into account, it will appear that the force 
on a stationary magnetic dipole due to a moving elec- 
tric monopole is not equal and opposite the force on the 
electric monopole due to the magnetic dipole. Newton's 
law of action and reaction will thus be violated, if hidden 
momentum is disregarded. 

There is also support in the literature [33|, HH for 
there being a hidden momentum associated with the 
intrinsically-magnetic electron, in the presence of an elec- 
tric field external to the electron. In any case, a quasiclas- 
sical analysis of particles with intrinsic magnetic moment 
will find Newton's law of action and reaction violated if 
hidden momentum is not included. The problem of vi- 
olation of the law of action and reaction led to the dis- 
covery of hidden momentum, and applies in the case of 
particles with intrinsic magnetic moments, when treated 
classically, just as it does in the case of classical current- 
carrying rigid bodies. Violation of the law of action and 
reaction would seem to lead inevitably to impossibility of 
constructing a consistent dynamical description, and so 
cannot be considered physical. Therefore for the present 
work it will be assumed that for a particle with intrinsic 
magnetic dipole moment /i, the total mechanical momen- 
tum must include a hidden momentum /x x E/c. The 
electron mechanical momentum is then 

-Pmcch,e = Pe + Phidden.e = P e + V e X E(r e )/c. (15) 

The rate of change of the electron mechanical momen- 
tum is accordingly 

-Pmech.e = P e + A* e x E(r e )/c + p, e x E(r e )/c. (16) 

Equating the rate of change of the mechanical momen- 
tum to the applied force obtains the rate of change of the 
electron kinetic momentum as 

P e = F e -fi e x E(r e )/c -fi e x E(r e )/c, (17) 
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where F e is the force applied to the electron. The 
latter will consist of the Lorentz force on the electron 
charge as well as forces due to field gradients acting on 
the electron intrinsic magnetic dipole moment, and on 
the electric dipole moment it aquires due to translational 
motion. 



C. Approximate Equation of Translational Motion 

The objective of the present analysis, as stated above, 
is to obtain an equation of motion of the total angu- 
lar momentum of quasiclassical Positronium atom model 
that is accurate to order /3 2 . The angular momentum 
contribution of the time rate of change of the kinetic 
momentum as given by Eq. (|17p is simply of the form 
r e x P e ; that is, each term in Eq. (fTTj) will contribute pro- 
portionately to the orbital angular momentum. There- 
fore, if any terms on the right hand side of Eq. (fT7|) are 
smaller than 1 times the largest term, their correspond- 
ing terms in the equation of motion of the orbital angular 
momentum may be neglected. 

In order to assess the relative magnitudes of the three 
terms on the right hand side of Eq. ([TT]) contributing to 
the time rate of change of the kinetic momentum, for the 
electron and positron in the quasiclassical Positronium 
atom, it will be useful to approximate the electric field 
as simply the instantaneous Coulomb field of the charge. 
The degree to which this is a valid approximation is con- 
sidered in detail in Appendix B. In particular, the electric 
field contributions due to the translational motion of the 
electron intrinsic magnetic dipole, and its relativistically- 
acquired electric dipole moment are not negligible to the 
present analysis and will be taken into account. However 
a small influence on a term that is already small may be 
negligble. 

It is shown in Appendix B that the electric field at the 
electron due to the positron, in the quasiclassical Positro- 
nium model with circular orbit, accurate to order of f3 2 , 
is 



E ( r e) = 72" + ^3 [ 3n ( V P X Mp) ' 71 - 2 ( V P X '( 18 ) 

(This neglects the acceleration field, which is an order 
/3 3 correction, and deviation between the instantaneous 
Coulomb field and the exact delayed Lienard-Wiechert 
velocity electric field, which is an order /3 4 correction. 
These errors are evaluated in the appendix.) 

Preparatory to considering what terms are relevant in 
the equation of motion of the kinetic momentum with 
an approximation for the field, it must also be consid- 
ered what errors may be introduced into the mechanical 
momentum itself (as given by Eq. (Tl3|) ) by any field 
approximation. This is essential because even though 
these errors are small in magnitude, they may be rapidly 
varying and so may contribute disproportionately to the 
equation of motion of the kinetic momentum. However, 



since some of the terms are already varying at the orbital 
frequency, any terms of smaller magnitude will not be- 
come larger after time differentiation if they do not vary 
faster than the orbital frequency. It is therefore assumed 
that there is no time variation faster than at the orbital 
rate in smaller-magnitude electric field terms in Eq. (|15[) . 

In order to determine what terms of the electric field 
and electric field time derivative should be retained when 
calculating the time variation of the kinetic momentum 
according to Eq. (TT71) . it is useful to assess the relative 
magnitudes of the terms under the assumption that the 
electric field is simply the Coulomb field. Since it is es- 
tablished that corrections to the Coulomb field and its 
first time derivative arise only at the order /3 2 , these cor- 
rections may be negligible in an equation of motion of 
the momentum where contributions of terms involving 
the hidden momentum are small compared to the kinetic 
momentum. Therefore Eq. (|17[) will be evaluated ini- 
tially under the assumption of a purely Coulomb electric 
field, although later a more exact expression for the elec- 
tric field will be used to determine the motion of the 
electron and positron orbits. 

The first term on the right hand side of Eq. (IT71) is the 
translational force on the electron. This force consists 
of the Lorentz force on the electron charge, plus field- 
gradient forces acting on the intrinsically-magnetic elec- 
tron. All non-Coulomb forces are evaluated in Section VI 
and shown to be of order /3 2 or smaller than the Coulomb 
force, in the Bohr model of Positronium ground state. 
Therefore similar considerations as support approximat- 
ing the electric field as the Coulomb field support ap- 
proximating the translational force as the Coulomb elec- 
trostatic force, in determining the approximate relative 
magnitude of the terms on the right hand side of Eq. 

(ED- 

The Coulomb field of a particle of charge, q, is 

^Coulomb = -3 • (19) 

With the electric field approximated as the Coulomb 
field, and the translational force approximated as the 
electrostatic Coulomb force, the first term on the right 
hand side of Eq. (TTT)) is then approximated as 

2 

F e ^eE^- E -^-. (20) 

The magnitude of the force on the electron is then 
given to the same order of approximation as 

\F e \* e l. (21) 

At the Positronium Bohr radius re = 2h 2 /me 2 this 
becomes 

4 2 6 2 

. _ , 9 e m e"m . . 

l^l* e2 ^r = W (22) 
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The magnitude upper limit of the second term on the 
right hand side of Eq. (|17p , assuming a Coulomb electric 
field and a circular orbit, is 



A* e x E/c 



e 3 h 



cr 3 y/2m 3 / 2 c 2 r 7 / 2 ' 



(23) 



where the Bohr magneton formula /ie = eh/2m 
and the result from Appendix A (Equation (|A2|| ) that 
the electron-positron relative velocity magni tude i n the 
circular-orbit Rutherford model is v — ey2/mr have 
been used. At the Bohr radius this becomes 



Me X E/C 



< 



16 



e 10 m 2 
c 2 W> 



(24) 



The ratio of the second term on the right hand side of 
Eq. (|17p magnitude to the translational force magnitude 
is thus bounded for the Bohr model ground state as 



Me X E/C 



< 



(25) 



where a = e 2 /he is the fine structure constant. This 
shows that the second term is nominally an order 1 
order correction to the change in mechanical momentum 
due to the Coulomb force alone. Therefore corrections to 
this term that are of order j3 or smaller may be neglected 
in the equation of motion of the mechanical momentum 
that is accurate to order 1 . 

The magnitude upper limit of the third term on the 
right hand side of Eq. (|17l) , assuming a Coulomb electric 
field and a circular orbit, is 



i . ^ ; i e/i e es 
Me x E c <-^ = — — 

cr* cr cm 

At the Bohr radius this becomes 



I Me x E/c\ < 



2h 2 



(26) 



An upper bound for the magnitude of the time rate of 
change of the spin vector, for the quasiclassical positro- 
nium atom in the Bohr ground state, is developed in 
Appendix D (see Eq. (|D3|) ). as 



The ratio of this bound to the first term magnitude is 



\fi e xE/c\ ^ 7 e 14 m 2 fcm- 



128 h s c 4 V h 4 



or 



lM e x E/c\ < 7 



7 



128 h 4 c 4 128 



(30) 



(31) 



The third term on the right hand side of Eq. (|17l) . 
assuming a Coulomb electric field and a circular orbit, 
is thus at most a factor of /3 4 as large as the first term 
that is the force on the electron. Therefore it may be 
considered negligible in the present analysis. 

In Appendix C it is shown that non- Coulomb terms in 
the electric field contribute to the rate of change of the 
electric field only at the order of f3 2 or below compared 
to the time rate of change of the Coulomb electric field, 
under the assumptions of the electron intrinsic magnetic 
moment equivalent to that of a classical current loop. 
Since the electric field time-derivative terms assuming 
the Coulomb field contributed to Eq. p7|) only at order 
f3 2 or below to the time rate of change of the kinetic 
momentum, non-Coulomb field contributions will not be 
needed in constructing an approximation for Eq. (|17p 
that is accurate to order f3 2 . 

Neglecting the third term on the right hand side of Eq. 
(TTT|) . and approximating the electric field in the second 
term as the Coulomb field, based on the considerations 
as described above, results in the following equation of 
translational motion for the dipole-carrying electron in 
the field of the co-circular-orbiting positron: 



cr 3 



M e x r, 



(32) 



where F e is the total translational force on the elec- 
tron. 

It is important to recognize that although it was as- 
sumed that the electric field was the Coulomb field in 
calculating the relative magnitude of the three terms on 
the right hand side of Eq. (fl~7|) . when the dynamical 
behavior is investigated the force on the electron charge 
will be based on the electric field approximation given 
by Eq. (fTSJ), which includes non- Coulomb terms due to 
the positron intrinsic dipole moment, and is accurate to 
order /3 2 . Forces on the electron due to its dipolar nature 
and the presence of field gradients are significant to the 
analysis and will be included as well. 



< 



so 



7 e 8 m 

l&Wd 2 ' 



i . t—i / i ( e 6 ™ \ 7 e 8 m 



7 



128 ft 8 c 4 



(28) 



(29) 



D. Motion of the Kinetic Orbital Angular 
Momentum 



The equation of motion of the electron kinetic orbital 
angular momentum around the center of mass based on 
Equation (|3"2"j) for the electron translational motion is 
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L e = r e x mti e = r e x F e ^r e x (p x r). (33) 

For the circular orbit and using the vector identity a x 
(b x c) = (c- a)b — (6 • a)c, and with v = r = 2f e , and 
noting that for the circular orbit -u • r = 0, this can be 
rewritten in what will turn out to be a more convenient 
form as 

L e= T e + ■ (l e )v, (34) 

where T e = r e x F e . For the positron it is similarly 
obtained that 

g 

The equation of motion of the total kinetic orbital an- 
gular momentum, L = L e + L p , is then 

6 

L = L e + L p = r e + t p + 7^—3 ' (Me ~ Mp))«- ( 36 ) 

In terms of spin vectors rather than intrinsic magnetic 
moments, with electron and positron g-factors g e — g p — 
2, and in terms of the total torque, 

e 2 

L = t-- — aj(r-(«e + «p))«, (37) 

where r = r e + r p is the total of the separate torques 
on the electron and positron orbits. 

With the general equation of motion of the kinetic or- 
bital angular momentum in hand in Eq. ()37|) . it remains 
to determine the total torque r, in order to evaluate the 
contribution of the kinetic orbital angular momentum to 
the total orbital angular moment, and in turn to the total 
angular momentum. 



1) Coulomb force. 

2) Biot-Savart force on the charges transiting the mag- 
netic field (consisting of the intrinsic magnetic field of 
the other particle and due to the motion of the other 
particle's charge). The former force will be referred to 
as i^Biot-Savart l- The latter force will be referred to as 

■^Biot-Savart 2- 

3) Force on the intrinsic magnetic moments due 
to anisotropy of the magnetic field (the Stern- Gerlach 
force). The Stern- Gerlach force due to anisotropy of the 
magnetic field due to the other charge motion will be re- 
ferred to as -Fstern-Geriach l- The Stern- Gerlach force due 
to anisotropy of the other particle's intrinsic magnetic 
field will be referred to as F S tcrn-Goriach 2- 

4) Force on the charges due to electric field induced by 
motion of the other particle magnetic moment. The force 
due to the electric field induced by translational motion 
of the other particle intrinsic magnetic moment will be 
referred to as Fvxfi- This force will also include the 
influence of the electric field due to the motion-acquired 
electric dipole moment of the other particle. (The influ- 
ence of the electric field induced by orientational motion 
of the other particle intrinsic magnetic moment will be 
shown to be insignificant to the present analysis). 

5) Force on the motion-acquired electric dipole mo- 
ments due to anisotropy of the electric field of the 
other charge and the other particle's motion-acquired 
electric dipole moment (the latter is shown negligible). 
The non-negligible former force will be referred to as 

-^Stern-Gerlach like- 

It is important to recognize that the forces act to 
change the total mechanical momentum of the particle 
that includes both the momentum of the moving mass 
and the "hidden" mechanical momentum of the intrin- 
sic magnetic moment in the electric field present at each 
particle's location in the Positronium atom. 



V. INVENTORY OF TRANSLATIONAL 

FORCES ON THE PARTICLES AND 
RESULTANT TORQUES ON THE ORBIT 

Preparatory to evaluating the rate of change of the ki- 
netic orbital angular momentum using Eq. (|37p , it will be 
useful to take inventory of all of the translational forces 
acting on the particles, as well as which of them will yield 
torques on the orbit. This will provide a road map for 
evaluating the total torque as needed according to Eq. 
(EU). 

A. Inventory of Translational Forces on the 
Particles 

The laboratory-frame translational forces on the par- 
ticles may be enumerated as follows: 



B. Summary of Torques on the Orbit 

The forces enumerated immediately above will lead to 
torques on the particle orbits if they have non-radial com- 
ponents with respect to the center of mass. 

The Coulomb force is exactly radial in the approxi- 
mation of the present analysis (which is shown in Ap- 
pendix B to be valid to order /3 3 ). Therefore it does not 
contribute any torque on either the electron or positron 
orbits. All of the other forces enumerated above may 
generate torques separately on the electron and positron 
orbits. It is desirable to identify them separately by par- 
ticle furthermore because they are added on the one hand 
in calculating the motion of the orbital angular momen- 
tum, and differenced on the other when calculating the 
motion of the orbital magnetic moment, as will be needed 
in a planned sequel. 
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VI. CALCULATION OF THE TORQUES ON 
THE ORBITS 



In this section the individual torques on the electron 
and positron orbits as identified above are explicitly cal- 
culated. (The reader not immediately interested in the 
details of the calculation of the torques may skip to Sec- 
tion VII, where the results of the torque calculations are 
summarized prior to their application in determining the 
equation of motion of the total orbital angular momen- 
tum.) 



A. Torque on Orbit Due to Biot-Savart Force in 
Magnetic Field due to Motion of the Other Charge 

The laboratory-frame magnetic field at the positron 
due to the motion of the electron may be written as 



B = ^(r e x(-r)). 



(38) 



where the minus sign of (— r) accounts for that r is 
defined as the displacement from the positron to the elec- 
tron, rather than towards the positron as the field point 
as needed here. 

The Biot-Savart force on the positron traversing the 
magnetic field is 



e 3n (to • u e ) - fi e 

Biot-Savart 2,p = ~V p X , [46) 



with n = —r/r here, and v p the proton velocity as 
measured in the electron rest frame. 

The torque on the positron orbit due to the electron 
intrinsic magnetic moment is then 



„ ' e 3n (to • fi e ) - Me i 

r = r p x F p = r p x ( -v p x ) . (44) 



The vector triple product involving u e can be ex- 
panded as 



r p x (v p x fi & ) = (r p ■ n e )v p - (r p ■ v p )fi e = (r p ■ n e )v p , 

(45) 

since the velocity and position vectors are orthogonal 
for the circular orbit so that the second term in the center 
vanishes. Similarly the vector triple product of Eq. (|44[) 
involving to yields that r p x (v p x to) = r p v p and so Eq. 
(|4"4")l becomes 



[3 (to • u e ) r p v p - (r p ■ n e )v p ] , (46) 



Biot-Savart l,p 



? >< (£ (*- >< ("*))) • ( 39 ) 



The torque on the positron orbit due to this force is 



^"Biot-Savart 1,; 



or 



(40) 



TBiot-Savart l,p = -T p X [f e (f p • r) - T (f p ■ T e )] 

° ^41) 

Since r p = —r/2, the cross product of r p with r van- 
ishes, leaving 



''"Biot-Savart 1,: 



T,x[f e ( V r)](|£). (42) 



This torque vanishes for the circular orbit, due to per- 
pendicularity of the radius and the velocity vectors, as 
does the similar torque on the electron orbit, and so need 
not be considered further herein. 



B. Torque on Orbit due to Biot-Savart Force on 
Charge Motion through Intrinsic Magnetic Fields 



''"Biot-Savart 2. 



P ~ S 



[2(r p ■ fx e )v p ] 



(47) 



Recognizing that r = — 2r p and v = —2v p , and re- 
placing the electron intrinsic magnetic moment with its 
equivalent in terms of intrinsic spin, obtains the final 
form for the torque on the positron orbit due to the Biot- 
Savart force in the electron intrinsic magnetic field as 



''"Biot-Savart 2,; 



2mc 2 r 3 



v [r ■ s e ] 



(48) 



The torque on the electron orbit due to the Biot-Savart 
force in the positron intrinsic magnetic field is found sim- 
ilarly to be 



''"Biot-Savart 2,- 



2mc 2 r 3 



v [r ■ s p ] . 



(49) 



Since the other type of Biot-Savart force-derived 
torque vanishes, the total torque on the orbit due to Biot- 
Savart forces is thus 



The force on the positron charge due to traversing the 
electron intrinsic magnetic field is 



''"Biot-Savart 



2mc 2 r 3 



v [r ■ (s e + s p )} . (50) 
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C. Torque on Orbit Due to Stern-Gerlach Force in 
Magnetic Field due to Motion of Positron Charge 

The force on the electron intrinsic magnetic moment 
due to the motion of the positron charge is given (in the 
laboratory frame) by 



Tstcrn-Gorlach l,p = T p X [V(/LX p ■ B)} = (n ■ fl p )v . 

(56) 



With r„ 



-r e , the total torque on the orbit due to 



the Stern-Gerlach 1 force is thus 



V (M e ' B ) = ^Stern-Gerlach l,e = V (-^ (/X e ■ (v p X r)fj , 

(51) 

which evaluates to 



-^Stern-Gerlach l,e = ^[/J e Xt) p - 3/X e • (t> p X Tl) «] . 

(52) 

In terms of the electron-positron relative velocity, v = 
—2v p , this becomes 



-Gcrlach l,e 



2cr 3 



[M e x v - 3/x e ■ (v x n) n] . 

(53) 



The torque on the electron orbit due to the Stern- 
Gerlach force on the electron intrinsic magnetic moment, 
in the magnetic field due to the motion of positron charge 
is then 



TStorn-Gorlach l,e = r e X [V(/X e • B)} = T e 



q P 

~2cr 3 



[H e x v] 
(54)' 

The vector triple product can be expanded as r e x(vx 
Me) = ( r e ' Ve) v ~ ( r e ' ^)M e = r(n ■ H e )v/2, since n = 
r/r = 2r e /r, and recognizing that the proton velocity 
and position vectors are orthogonal for the circular orbit 
so that the second term in the center vanishes. It is thus 
obtained that 



^"Stcrn-Gcrlach l,e 



q P 

Acr 2 



(n- n e )v. 



(55) 



The similar torque on the positron orbit is 



''"Stern-Gerlach 1 



4cr 2 



(n-(p e -n p ))v. 



(57) 



This is rewritten in terms of intrinsic spins rather than 
intrinsic magnetic moments, and replacing n with r/r, 
as 



TStorn-Gorlach 1 = ~ 4mc 2 r 3 ( r ' ( Se + S p)) V ' ( 58 ) 

D. Torque on Orbit Due to Stern-Gerlach Force 
due to Intrinsic Magnetic Field of Other Particle 



The force on the electron due to anisotropy of the mag- 
netic field due to the positron intrinsic magnetic moment 
is given by 



-^Stern-Gerlach 2,e 



(59) 



where 



B ^ P - 



3n (n • n p ) - n p 



(60) 



The force on the electron is then 



^Stern-Gerlach 2,e — V ^Me 

which evaluates to 



3n {n- n p ) -Mp" 



(61) 



^Stern-Gerlach 2,e = - (j^j (/* e ' Mp) - 3 (Me Mp) (^f) + 

(^) (^e ' n) X - [n p - n{n p ■ n)] + (J^j (n ■ p p ) J [/x e - n(p e ■ n)] . (62) 
The torque on the electron orbit in the laboratory frame is then given by 

Tstcrn-Gcrlach 2,e =T e X [V(Me ' B)] =T e X [fl p (fl e ■ Tl) + M e (" " Mp)] ■ (63) 

Similarly, for B at the positron due to the electron intrinsic field, 
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T"Stc 



-Gcriach 2, P = r p x [V(/x p ■ B)] = r p x [fi e (fi p ■ (-n)) + ((-n) • /ij Mp] f |rj 
or, since for equal- mass particles r p = — r e , 



(64) 



TStcrn-Gorlach 2,p = T p X [V(/ip ' #)] = X [/^ (/X p • Tl) + (n ■ fJ, e ) /X p ] f -j 

The total torqe on the orbit due to this category of force is then (with r = r e — r p = 2r e ) 



TStern-Gcrlach 2 = V X [fl p (fl e ■ Tl) + (Vl ■ jU p ) /ij , 



or, in terms of the intrinsic spins, 



(65) 



(66) 



TStern-Gcrlach 2 = ~T X [s p (s e • Tl) + (tl ■ S p ) S e ] 

This can be put into a form that can be more directly used later as 



3e 2 



(67) 



^"Stcrn-Gcrlach 2 



— [s p x 3n (s e • n) + s e x 3n (n • s p )} 
m err 



(68) 



E. Torque on Orbit Due to Electric Force due to 
Motions of Intrinsic Magnetic Dipoles 

A moving magnetic dipole gives rise to electric field due 
to the time- varying magnetic field it causes. The electric 
field due to the motion of each magnetic dipole-carrying 
particle acts on the other particle's charge. Also, a trans- 
lating magnetic dipole must acquire an electric dipole 
moment [35| . This also gives rise to electrical forces act- 
ing on the other-particle charge. These forces give rise to 
orbital torques if they have components that are nonra- 
dial around the center of mass. 

The total electric field due to translational motion of 
the magnetic dipole, for the Positronium atom model 
Bohr ground state, and accurate to within a factor of 
/3 2 times the Coulomb electric field strength, is found in 
Appendix B4 to be 



J Vx/jb 



■ (S((v x fi) ■ n)n - 2(v x n)) . (69) 



The force accurate to order (v/c) 2 acting on the elec- 
tron due to the circular orbiting intrinsically-magnetic 
positron is then 



Fvxfi = -— (3((u p x fi p ) ■ n)n - 2(v p x fj, p )) . (70) 
The torque on the electron orbit is thus 



T Vx(Jb,e — r e 



— g2(Vp X fl p ) 



(71) 



or, returning to the convention that v is the electron 
velocity relative to the positron, and in terms of the 
positron spin rather than the intrinsic moment, 



T~VxfJ,,e 



r x 



-(v X s p ) 



(72) 



2mc 2 r 3 

The torque on the positron orbit is found similarly to 



be 



TVxfl,e = -T 

and the total torque, 



2mcr 3 



(v X s e ) 



(73) 



[v x (s e + s p )} . (74) 



Vx L l 2mcr 3 
For the circular orbit this can be simplified as 

e 2 
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F. Torque on Orbit due to Stern-Gerlach-like Force 
on Electric Dipole 

A magnetic dipole m translating with velocity v ac- 
quires an electric dipole moment, d, as [35j 



d = v x m/c. 



(76) 



The intrinsically-magnetic electron thus acquires an 
electric dipole moment d e due to its translational mo- 
tion given by d e = v e x fl e /c. 

The force on an electric dipole d in a electric field E is 



F = V{d-E) = (d-V)E, 



(77) 



where the second equality applies provided that Vx£ 
vanishes or is negligible as in the present application. 

For the Coulomb electric field due to the positron 
charge, acting on the electron, this becomes 



(d e -V 



(78) 



With 



d 3(d ■ n)n 



the force on the electron becomes 



F e = [(v e x n e ) - 3((v e x fi e ) ■ r)r] . 
The torque on the orbit due to this force is thus 



T e = r x F e = ^r e x [{y e x /xj] 



(79) 



(80) 



(81) 



or, in terms of the spin instead of the magnetic mo- 
ment, and with q p = —q e = e, and r = 2r e and v = 2v e , 



t p = -- 



T r x [v x s e ] . 



4mc 2 r 3 

The torque on the positron orbit is similarly 

r x [(v x s p )] . 



(82) 



4mc 2 r 3 



(83) 



The total torque on the orbit due to this type of force 
is thus 



T rx[vx (s e + s p )} . (84) 



4mc 2 r 3 

For the circular orbit this torque can be rewritten as 



^"Stcrn-Gcrlach-likc 



-4^3 V t r ' ( S - + S f)] • ( 85 ) 

G. Negligibility of Electric Dipole-to-Dipole Force 
and Torque 



A Stern-Gerlach-like force and corresponding torque 
exists due to anisotropy of the electric field of the motion- 
acquired electric dipole moment of one particle, acting on 
the motion-acquired electric dipole moment of the other. 
This torque is shown to be negligible to the present anal- 
ysis. 

Using the first equality of Eq. (|TT[) obtains for the 
electric field due to the positron electric dipole moment 
acting on the electron the force 



F = V [d, 



3(n • d p )n — d f 



(86) 



This force is formally similar to the Stern-Gerlach force 
on the electron intrinsic magnetic moment in the intrinsic 
magnetic field of the positron (as given by Eq. (I5T1) ). It is 
weaker by a /3 2 factor however due to the v/c in Eq. ([76]) 
for the acquired electric dipole moment and is therefore 
negligible in the present analysis. 



VII. MOTION OF THE TOTAL ORBITAL 
ANGULAR MOMENTUM 

The total orbital angular momentum, referred to here 
as Ltotai, consists of the total kinetic orbital angular mo- 
mentum, that is here represented as L = L e + L p , the 
total hidden mechanical orbital angular momentum, rep- 
resented as .^hidden, and the total electromagnetic field 
orbital angular momentum. The latter is due to the elec- 
tron and positron intrinsic magnetic fields crossed onto 
the electric field of the other particle, and will be rep- 
resented as ificid- The equation of motion of the total 
orbital angular momentum is then expressible as 



J total 



'hiddc 



afield- 



(87) 
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A. Rate of Change of the Kinetic Orbital Angular Momentum 

In this section the rate of change of the total kinetic orbital angular momentum, that is, the total orbital angular 
momentum of the translational motion of the electron and positron masses, is calculated. 
It was determined previously (Equation (j3"7| ) that 

e 2 

L = L e + L p = t - (r • (s e + s p ))v 

where r = r e + t p is the sum of the torques on the electron and positron orbits. From the results of the previous 
section, the sum of the non-negligible torques is 

T = TBiot-Savart + T Stern-Gcrlach 1 + Tstcrn-Gcrlach 2 + Ty X/J, + T Stern-Gerlach like 

where, from Equations O, ([55]). ([75]). (|55|) , and with n = r/r, 

e 2 

TBiot-Savart = ~ 2m C 2 r 2 " ^ ' ^ + ^ 



T Storn-Gcrlach 1 — ~ "J 5-5 {n ■ (s e + S p )) V 

4mrr 



Tstorn-Gcrlach 2 = 2 2 5 [Sp X 3n (s e ■ Ti) + S e X 3ri(n • Sp)] 



TVX/X "~W3 [(»" ' («e + «p))«] 



TStern-Gerlach like — — ~. 9 9 ^ [« 1 («e + s p)] 

Evaluating the sum of the torques according to Eq. ([88]) obtains the total torque as 

3e 2 e 2 

t = - ~ 5-5- (" • (s e + Sp)) v + „\ [sp x 3n (s e • n) + s e x 3n (n • s p )l . (89) 

Zmc z r z r''m z c 

The rate of change of the total kinetic orbital angular momentum evaluated according to Eq. (|3"T1) is thus 

2e 2 e 2 

£ = 5-^ (« ■ (s e + s p )) v + — 5— [s p x 3n (s e ■ n) + s e x 3n (n • s„)] . (90) 

mc r r a m z c 

It will turn out to be useful to understand the relative magnitude of the two terms summed on the right hand side of 
Eq. (IM1) . as well as the overall magnitude of L. This will aid in demonstrating the negligibility of some contributions 
to the motion of the total angular momentum, and so simplify the analysis. The first term on the right hand side of 
Eq. ([90]) is bounded as 



2e 



2 



(n ■ (s e + s p ))v 



< 2e2 n ^ = 2 ^ h , (91) 



The second term on the right hand side of Eq. is bounded as 
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[s p x 3n (s e ■ n) + s e x 3n (n ■ s p ) 



< 



6e 2 



r 3 m 2 c 2 4 



(92) 



The ratio of the second term bound to the first is then 



6e 2 



2y/2e 3 h 



r 3 m 2 c 2 4 I c 2 TO 3/2 r 5/2 / 4\/2 \em 1 / 2 r 1 / 2 
At the Positronium Bohr radius re = 2h 2 /(e 2 m), this becomes 



(93) 



4^ \ern 1/2 



2ft 2 



-1/2 



4^ Vem 1 / 2 / I y2ft 



,1/2 



(94) 



The two terms on the right hand side of Eq. (|90|) being of similar magnitude are thus both significant to the rate of 
change of the kinetic orbital angular momentum, for the quasiclassical Positronium atom in the Bohr ground state. 



B. Rate of Change of Hidden Orbital Angular Momentum 

With the hidden momentum of a magnetic dipole of moment m as given previously asmx£/c, the hidden angular 
momentum around the center of mass, of the electron intrinsic magnetic moment in the Coulomb field of the positron, 
is defined as 

1 e 

£hidden,e = ~T e X (fl e X E) = —»T e X X r) . (95) 

C CT' 3 

(Since the hidden orbital angular momentum magnitude is about a factor of /3 2 smaller than the kinetic orbital 
angular momentum magnitude, there is no need to include higher order terms in the description of the electric field 
here.) Expanding the vector triple product and with r = 2r e obtains that 



J hiddcn.i 



2cr 3 



[r 2 ii e -r(n e -r)]. (96) 



For the circular orbit (where r = 0), 



rf£hiddcne,e _ l h . e = [ r 2^ — v (fi ■ r) — r ((/it e ■ r) + (p, e ■ »))] . (97) 
at 2cr A 

Rearranging and in terms of the electron spin instead of the intrinsic magnetic moment, 

e 2 e 2 

ihiddcn, e = i v ( s e -r)+r(s e - v)] - 5- [s e -n(s e - n)] . (98) 

Zmc r° Zmc r 

Now, the second term in the brackets on the right hand side can be further evaluated by substitution for s e as given 
by Eq. @ , but it is worthwhile instead to consider the relative magnitude of the two terms being differenced in Eq. 
(|98l) . As has already been used in Section IVC, an upper bound for the magnitude of the time rate of change of the 
spin vector, for the quasiclassical positronium atom in the Bohr ground state, is developed in Eq. (|D3[) in Appendix 
D. Taking r to be the Bohr radius then obtains that 



2mc 2 i 



[b e - n (b e ■ n) 



< 



e 2 e 2 m 7 e s m 



7 em 



2mc 2 H 2 32 h 4 c 2 64 h 6 c 4 



(99) 



An upper bound on the first term on the right hand side of Eq. (|9"5|) can be developed as 



19 



[v (s e ■ r) +r (s e ■ v)} 



< 



e% /2 h 



5 V2 



2mc 2 r 2 2mc 2 r 2 m 1 / 2 r 1 / 2 2 c 2 TO 3/2 r 5/2 4 ' 



which at the Bohr radius becomes 



(100) 



[v (s e ■ r) + r (s e ■ v)] 



< 



e 3 \/2 h ( 2h 2 



2mc 2 r 3 

The relative magnitude of the two bounds can now be evaluated as 



-5/2 



I6c 2 h 4 ' 



(101) 



e 8 m 



7 e 12 m 



64 H 6 c 4 V I6c 2 h 4 



7 e 12 m {I6c 2 h 4 



64 h 6 c 4 



7 e 4 



4h 2 c 2 



(102) 



The time rate of change of the electron hidden orbital angular momentum can now be written accurate to one part 
in 1//3 2 as 



■'hidde 



2mc 2 r 3 



[r (s e ■v) + v(s e - r)] 



(103) 



The rate of change of the hidden orbital angular momentum of the positron intrinsic magnetic moment in the 
electric held of the electron is evaluated similarly to obtain 



iddcn,£ 



2mc 2 r 3 



[r (s p -v)+v(s p - r)] 



(104) 



The rate of change of the total hidden orbital angular momentum is then (accurate to one part in 1//3 2 ) 



ddc 



2mc 2 r 3 



v((s e 



r)] 



(105) 



In order to confirm that Eq. (|105l) is sufficiently accurate for the present purpose, the relative magnitude of the 
rate of change of the hidden orbital angular momentum to that of the kinetic orbital angular momentum as given by 
Eq. (|90[) can be evaluated. Since the two terms on the right hand side of Eq. (|9"0|) were found to be of the same order 
of magnitude (at the Bohr ground state), either can be taken as representative of the general magnitude of the rate 
of change of the kinetic orbital angular momentum for the Bohr ground state. The ratio of the bound on the larger 
term of the hidden orbital angular momentum motion, from Eq. (|100[) . to the bound on the first term on the right 
hand side of Eq. (|9"0|) , as given by Eq. (I^Tt , is 



/ 2y/2e 3 h \ _ 1 



! V2 



(106) 



This result establishes that the magnitude of the rate of change of the hidden orbital angular momentum is similar 
to that of the kinetic orbital angular momentum. 

It is also of interest to compare the magnitude of the term involving s e , in the motion of the electron hidden orbital 
angular momentum of Eq. (|98l) , to the magnitude of the contribution to the motion of the total angular momentum 
attributable to the Thomas precession. The contribution of the Thomas precession to the motion of the electron spin 
vector in the quasiclassical Positronium atom as given by Eq. (j9]), evaluated at the Bohr radius, is bounded as 



< 



2h 2 



32c 2 h 4 ' 



(107) 



The ratio of the magnitude of the term involving s e in Eq. (|98[) to the magnitude of the contribution to the motion 
of the total angular momentum attributable to the Thomas precession is then 
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7 ' j l» 2 - (108) 



32c 2 fr 4 j 2 H 2 c 2 

Since (3 is the same order as a here, neglecting the term involving s e in Eq. for the hidden orbital angular 

momentum rate of change will introduce only order /3 2 errors compared to the effect of the Thomas precession, in the 
rate of change of the total orbital angular momentum. 



C. Rate of Change of Field Angular Momentum 



The total angular momentum can be expected to be a conserved quantity only if the field angular momentum is 
included [36| . For a system consisting of a stationary current distribution and a particle with electric charge q, and 
with r the displacement of the charged particle from the curent distribution, there is field angular momentum around 
a fixed point r\ evaluated as [37j 



ifieid = —J{r- ri )x(ExB) d 3 r = (q/c)(r - r x ) x A(r), (109) 

where A(r) is the vector potential of the current distribution evaluated at the location of the charged particle. 
If the vector potential in Eq. (|109p is taken to be that of the positron magnetic dipole moment, the vector potential 
at the location of the electron is 




where r in this case being the displacement of the electron charge from the positron is in agreement with the 
convention of the present analysis. 

With the electron the charged particle, q = —e, and taking the point T\ to be the center of mass then obtains that 
r — ri = r e (since r = r e — r p and the displacement of the center of mass from the positron position is — r p ), and so 
the field angular momentum according to Eq. (|109[) becomes 



***** = ~^3 r - X (»P X r ) = -^3 r X (»P Xr )= V 2 »P ~ r (*V r )] > ( U1 ) 

where the first equality is based on that r — 2r e . Noting that the expression of Eq. (|111[) for the total field angular 
momentum of the positron charge Coulomb field crossed with the intrinsic magnetic field of the electron is formally 
identical to the expression of Eq. for the hidden orbital angular momentum of the electron, allows the rate of 
change of the field angular momentum due to this contribution accurate to one part in 1//3 to be written immediately, 
based on Eq. (|103[) . as 



£fieid, P = 7; 5-^7 [r {s p -v)+v(s p - r)l . (112) 

Strictly, since Eq. (|109|) assumed the current distribution was stationary, Eq. (|112|) is the rate of change of field 
angular momentum in the positron rest frame, which is rotating due to Thomas precession. However, applying the 
standard transformation |38[ for the rate of change of a vector from a rotating to a nonrotating frame, it is found 
that this changes the expression of Eq. (jlll|) only by about one part in 1//3 2 = 1/a 2 here and so is negligible to the 
present analysis. 

The rate of change of the total field angular momentum is found by summing the contribution of Eq. (|112l) and 
the similar quantity for the field angular momentum of the positron Coulomb field crossed onto the electron intrinsic 
magnetic field, with the result 



-Afield = 



2mc 2 r 3 



[r ((s e + s p ) ■ v) + v ((s e + s p ) ■ r)} . 



(113) 
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D. Rate of Change of the Total Orbital Angular Momentum 

The rate of change of the total orbital angular momentum is evaluated as 



From Eqs. (jTU5|) and (|TT5|) . 



iddcn T" -Afield 



•'total — L + ^hidden + -Afield- 



[n ((s e + S p ) ■ v) + V ((s e + Sp) ■ Tl)] 



Substitution of the results of Equations (|90| and (|115j) into Eq. (|114| obtains 

e 2 e 2 

i-totai = o 2 5 [s P x 3n (s e ■ n) + s e x 3n (n • s p )] H [r ((s e + s p ) • v) - v ((s e + s p ) ■ r)} 

m err mc r 

With the vector identity ax (b x c) = (a ■ c)b — (a ■ b)c, 



itotai = 9 9 . [sp x 3n (s e • n) + s e x 3n (n • s p )] H ^ [(s e + s P ) x (r x v)] 



(114) 



(115) 



(116) 



(117) 



or (note r x mt) = 2L here), 



^total 



[s p x 3n (s e • n) + s e x 3n (n • s p )] 



2e 2 



[(s e + s p ) x L] . 



(118) 



VIII. MOTION OF THE TOTAL ANGULAR MOMENTUM 

With total angular momentum J = s e + s p + -Ltotai an d J = s e + s p + L tota \, it is obtained from Eqs (fTU|). and 
(jnSJ) that 



2 9 S 

m z c z r A 



3n (n • s p ) — s p 



2 9 3 P 



3n(n • s e ) - s e 



[sp x 3n (s e ■ n) + s e x 3n (n ■ s p )] + 



2e 2 



[(s e + Sp) x L\ , 



(119) 



which reduces to 



- [(8 e + s p ) x L] 



(120) 



This implies that the total angular momentum is a constant of the motion if s e j_ = —s p j_, where s e ±, s p ± are the 
components of the electron and positron spin angular momenta perpendicular to the total orbital angular momentum, 
but not otherwise, if the spins and orbital angular momenta are not all aligned. In general, the amount of nonconser- 
vation of angular momentum is exactly that attributable to the Thomas precession, as given by the sum of Eqs. (|11[) 
and ([12"]), 

The result of the present analysis, that total angular momentum is nonconserved due to Thomas precession, 
is surprising and perhaps distressing if correct. Therefore, it's worth supposing this result simply incorrect and 
considering what type of changes would be required to arrive at the expected result of J = 0. This may help in 
finding the expected errors. 

One simple change that will result in total angular momentum conservation is to change the 3/2 factor on L in the 
two spin equations of motion to 2. But, the value of 3/2 resulted as a collection of a factor of unity with a factor 
of one-half, the latter of which is the one-half that is recognized as the "Thomas factor" that resolved the spin-orbit 
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coupling anomaly of the anomalous Zeeman effect. This result is very well established. So, keeping the term with the 
one-half unchanged, the first term in the sum, that derives from the magnetic field at one particle due to the orbital 
motion of the other, would have to become 3/2 instead of unity. It's difficult to see how a "3" could creep into this 
term, however since the electron-positron displacement is simply r — 2r e — 2r pi and the L here has resulted from 
the particle velocities crossed onto the local Coulomb field. None of these involve any factors of 3. 

Alternatively, if the spin equations of motion derived here are taken as correct, then the other contributions to 
the total angular momentum must be in error, if Eq. (|120l) is incorrect. Since the field angular momentum does not 
vanish generally (if Eq. (|113[) is correct, or correct to within a factor), this would imply that the total mechanical 
momentum is not itself a constant of the motion (a prima facie interesting result, in itself, which is also true according 
to the present analysis). However, obtaining the expected result of J = here is not a matter of simply increasing or 
decreasing the field angular momentum or hidden orbital angular momentum by a factor or factors, or finding sign 
errors or missing factors on the individual torques acting on the orbit, alone. Rather, it would require at least two 
separate errors that when corrected arrived at the expected result. This is true because the form of the contribution 
of the Thomas precession, which as noted is the entire angular momentum nonconservation of Eq. (|120[) . involves the 
cross product of the spin vectors with the kinetic orbital angular momentum vector. Therefore it can be written as a 
vector triple product when the kinetic orbital angular momentum is expanded as 2L — r x mv, and then using the 
vector identity a x (b x c) = (c ■ a)b — (b ■ a)c obtains alternatively 



j = (]) [(s e + s p ) x (r x «)] - -Jr^ (]) [(« • (s e + Bp)) r - ((s e + s p ) -r)v]. (121) 

The form of the right hand side of this equation stands in contrast to the form of the equations for the hidden 
orbital or field angular momenta, Eqs. (|105l) and (111 3[) . involving a difference rather than a sum of otherwise similar 
terms. The minus sign here follows directly from an indisputable vector identity, while the plus sign in the equivalent 
position of the hidden orbital and field angular momenta follows from the also indisputable law of the differentiation 
of a product. If these signs are taken to be fixed then multiplying the hidden orbital or field angular momenta by 
any factors will not arrive at angular momentum conservation. It will have to be combined with a change in the 
total torque on the orbit, or the other term in the equation of motion of the kinetic orbital angular momentum of Eq. 
(|90| , that accounts for the presense of the hidden momentum. This is not impossible but a determined search by the 
author has turned up no such errors, and in fact it involved a very careful analysis simply to arrive at the situation 
reported here. 

By this time the reader is probably seeing the futility of this exercise, which is mentioned primarily to demonstrate 
that the possibility of an error or errors has been seriously considered and that errors have been sought for. Obviously 
the error search has been fruitless, or the result unhappy as it is would not be so presented. Therefore the reader 
is entreated to attempt to find the obstensible error, any identification of which will be welcomed. Nonconservation 
of the total angular momentum here is perhaps only a distraction from the equally interesting result that there is a 
complex dynamic involved in the electromagnetic two-body problem with spin that seems not to have been noticed 
previously. This dynamic would still be present even were the total angular momentum conserved, because the kinetic 
angular momenta and the motion of the spins would then have to be counter to the hidden and field angular momenta, 
which don't vanish independently. The total electric and magnetic dipole moments would thus not be constrained to 
be generally stationary, separately from the well-recognized dipole acceleration due to orbital motion of the separated 
charges, which clearly causes decay of the classical atom in the absence of spin. That there are additional dipole 
acceleration components not previously considered invalidates, at least until their full consideration, the conclusion 
that the classical atom must of necessity and unconditionally decay radiatively on a short time scale. 

Alternatively, considering that perhaps the result herein is correct, it can be noted that although the total angular 
momentum is nonstationary, the total magnetic moment of the quasiclassical positronium atom is not so constrained. 
This was noted previously in the quasiclassical Hydrogen atom [39j . where although the total angular momentum 
was not a constant of the motion, the total magnetic moment could nonetheless be stationary. The situation in the 
quasiclassical Positronium atom however is much richer, owing to the more rapid motions of the spins remarked about 
near the end of Section III. This motion may conceivably give rise to both magnetic and electric dipole radiation, 
since the particles acquire electric dipole moments due to their motions. In quasiclassical Hydrogen in the simple 
circular orbit, for example, the oscillation of the acquired electric dipole moment can cancel the oscillation of the 
dipole due to charge separation, at an orbital radius that is between the nuclear and atomic scale. Therefore such 
radiative effects and their reaction may be physically significant and potentially extend the applicability of classical 
physics further into the quantum domain. 
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IX. CONCLUDING REMARKS 



As described in the introduction, the 1927 analysis by L. H. Thomas found no angular momentum nonconservation 
in hydrogenic atoms due to the relativity precession. A detailed analysis showing how unawareness of the existence 
of hidden momentum can mask the angular momentum nonconservation and that this led to Thomas overlooking the 
nonconservation is provided in [ipf . The present analysis can thus be seen as in agreement with expectations based 
on the results of the analysis approach taken by Thomas, if performed according to modern textbooks recognizing 
the need for accounting for hidden momentum. 

It's also worth observing that although Thomas's 1927 analysis of hydrogenic atoms found no angular momentum 
nonconservation due directly to the "relativity precession," it did not either obtain strict angular momentum conser- 
vation, except in the trivial case of alignment of the spin and orbit angular momenta. Rather, in Thomas's analysis, 
only the "secular," i.e., orbit averaged, total angular momentum was a constant of the motion for general relative spin 
and orbit orientations. This is in contrast to the present analysis. Despite angular momentum nonconservation due 
to the Thomas precession in certain configurations, the present analysis obtains exact angular momentum constancy 
under the nontrivial conditions described, without need to average over an orbital period, but requires accounting for 
both the hidden and field orbital angular momenta, as well as the kinetic orbital angular momentum. 



Appendix A: Circular Orbit Rutherford and Bohr 
Models of Positronium 



Suppose a classical point-charge electron and positron, 
without intrinsic magnetic moments, are co-circular or- 
biting opposite their common center of mass. Equating 
the centripetal acceleration of electron circular motion 
with the acceleration due to the Coulomb attraction from 
the positron obtains 



so 



(Al) 



L e = L p = e^mr e ( - ] L. 



and since L 



L = e^/m e r e L. 



(A5) 



(A6) 



Alternatively, in terms of the electron-positron separa- 
tion, 



where m and v e are the electron mass and velocity 
measured in the center-of-mass frame, r e is the electron 
distance from the center of mass, and r is the electron- 
positron separation. The positron mass is also m. Also, 
r — r e +r p , where r p is the positron distance from the cen- 
ter of mass, and from the definition of the center of mass 
and since the electron and positron masses are equal, 
Te = r p , and r — 2r e — 2r p . Also, v p = —v e . If 
v is the electron velocity relative to the positron, then 
v = v e — v p , or v = 2v e . Through substitutions in Eq. 
(|A1[) it is then obtained that 



V2e 



(A2) 




L. (A7) 
Setting L to h obtains the Bohr radius of positronium 



as 



2h 2 



rB = 



1.05 x l(T s cm. 



(A8) 



(Note, this is the electron-positron separation, not 
strictly the orbital radius, but is the more convenient 
quantity for determining scale.) 

The electron velocity at the Bohr radius is, from (|A3|) 



The electron velocity as measured in the laboratory 
frame (assuming the Positronium atom center of mass is 
stationary there) is 



V2 



1 



2 2 y/mr 2 y/mr e 
The electron orbital angular momentum is 

L e = r e x mfl e = L e L = mr e v e L. 



(A3) 



(A4) 



V2 _e_ (e^m> 1/2 
2h 2 



and 



2 h ' 



(A9) 



(A10) 



where a = e 2 /(he) « 1/137 is the fine structure con- 
stant. 
Also 



24 



v e ps 1.10 x 10 8 cm/sec pa 0.004c, (All) 
which obtains 

/3 e 2 = (w e /c) 2 ps 1.33 x 10~ 5 , (A12) 

and 

7 - 1 ee (1 - f3 e 2 y 1/2 - 1 PS 6.66 x 10~ 6 . (A13) 



to the particle electric charges and due to their intrin- 
sic magnetic moments. An electron or positron in the 
Positronium atom will feel an electric field due to the 
charge of the other particle. There is also an electric field 
contribution due to the motion of the intrinsic magnetic 
moment, but this will be seen to be of magnitude or- 
der (v/c) 2 compared to the Coulomb field in the present 
analysis. As will also be shown, in the atomic scale the 
magnitude of the hidden momentum of an electron in the 
Coulomb electric field due to a proton is of order {v/c) 2 
compared to the electron kinetic momentum magnitude. 



Appendix B: Electrogmagnetic Field of 
Quasiclassical Positronium 

The electromagnetic field in the vicinity of the quasi- 
classical Positronium atom consists of contributions due 



The Lienard-Wiechert field expressions in three- vector 
notation are \4M 



E(r,t) 



n- (3 



7 2 (1 - (3 ■ nf R 2 



n x ( (to — (3) x /3 
{\-f3-nf R 



(Bl) 



and 



1. Electric Field due to a Circular-Orbiting Charge 



B{r,t) = [nxE] I 



(B2) 



where, if r$ is the displacement from the particle to 
a field point, then R = \ro\, to = tq/R. Also (3 = v/c, 
where v is the particle velocity. The subscript "ret" refers 
to that the quantity in the brackets is evaluated at the re- 
tarded time. The Lienard-Wiechert field expressions are 
the exact fields for a possibly moving point charge, but 
in the present application it must be taken into account 
that the particles have intrinsic magnetic fields as well. 
This will cause electric fields at least due to translational 
motion of the particles. However for the time being this 
will be disregarded as the electric field form needed due 
to the charges is determined. 



The first term on the right hand side of Eq. (IB1[) is 
called the velocity field. The second is called the accelera- 
tion field. In the present application, for the electric field 
experienced by one of the mutually-circularly-orbiting 
charges due to the other, the magnitude of the accelera- 
tion field is less than (v/c) 2 times the magnitude of the 
velocity field. It can also be seen by inspection that the 
magnitude of the field difference from the Coulomb field 
of the particle (that is, the electric field in the rest frame 
of the particle) due to motion is small when j3 « 1. 



Consider a positron and electron classical point charges 
in a classical circular orbit around their common center of 
mass. The positron and electron are assumed bound to- 
gether through the Coulomb force into a perfectly circu- 
lar Keplerian orbit, such that the particles are exactly op- 
posite each other in the orbit at all times to a laboratory- 
frame observer, and the particle velocity vectors are at all 
times perpendicular to the line between them. From the 
conditions of orbit circularity and direct opposition it is 
apparent that the particle separation r is constant, and 
so [r(t)]ret = r (t) = r - From the perpendicularity of the 
velocity to the line between the particles, it is apparent 
that (3 ■ n = at both the present and retarded times. 

It should be borne in mind that the perfectly circular 
orbit is only approximation, because forces other than the 
Coulomb attraction will cause deviations from circularity. 
It seems reasonable though that since the non- Coulomb 
forces are relatively small, the deviations they cause from 
perfect circularity will also be small, and can be assessed 
later as perturbations. 

The positron velocity magnitude in the circular 
Coulombic orbit is (see Eq. (IA3|) ) 



ey/2 

v p = o ; 

Zy/mr 

Choosing a coordinate frame to suit, suppose 



(B3) 
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that r is the electron-positron separation, and n is a unit 
r r vector in the direction of the electron from the location 

r pW = = 2 ( sinw ** ~ cosujt y) ( B4 ) of the positron. Then 

where the definitions have been retained from above 

I 

v p {t) = r p (t) = w{r/2) (cosuitx + sinwi) = v p (cosutx + sinwiy) = v p v(t), (B5) 

where 



v p e 



V2 



• p 

The velocity at the retarded time t — r/c = t — D is then 



_1 v ey/2 



r rri}-l 2 r 3 l 2 



(B6) 



KWlet = v p ( cos ( w ( f + sin (w (t - £>)) y) . (B7) 

Now 

(ui 2 D 2 \ ( uj 3 D 3 \ 
1 — J+smutlujD — J (B8) 

and 

sin (u> (t — D)) = s'mujt coswD — cos tut sin cjD w sinwt ( cjD — J — coswt I 1 — J . (B9) 

With ojD = (v/r)(r/c) — (3, and retaining terms to order /3 2 , 

cos(w(i- D)) w (1 - /3 2 )coswt + /3sinu;f = (1 - /3 2 ) coswi + /3coscu(t - n/{2u)) (BIO) 

and 

sin (w (£-£>)) w (l-/3 2 )sina;i-/?cosu;*; = (1 - /3 2 ) sinwi + /3smuj(t - 7r/(2w)), (Bll) 

so 

K W] rct ~ M*) + - r / 4 ) = + wD» P (i - 7r/(2w)), (B12) 

where T = 2-k/lo and 

[n(t)] iet « n(t) + ^«(t - 7r/(2w)). (B13) 
Alternatively (with sin(# — 7r/2) = — cos 6* and cos(8 — n/2) = s'md), 

[n(t)} Tet ^n(t)+wDv(t) (B14) 

and 

KWU « w (^/ 2 ) [*(*) - . (B15) 
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or 



(B16) 



For the circular orbit and evaluation of the field due to the positron at the position of the electron, v p (t)-n(t) = 0, so 
the Lienard-Wiechert electric field becomes at the electron (for the time being the subscript on the velocity indicating 
the positron will be omitted) 



E(r,t) 



n x ( (n — (3) x (3 



(B17) 



Note, the vector n here is from the source to the field point. This will be consistent with the definition above 
provided that the field is due to the positron and the field point is at the electron. The first term on the right hand 
side is the velocity field term and using Eqs. (|B14[) and (|B16I) evaluates to 



E ve i(r,t) = e 



n (3 



2 2 



n(t) - f3(t) 



2 2 



ecuD 



v(t) +ujDn(t) 



2 2 



(B18) 



or (to the order of the approximation) , with the Coulomb electric field due to the positron charge given as .Ecoui 
en/r 2 , and replacing (3 with v/c, 



E ve i(r, t) -E Cou \ 



eujrv{t) 
cy 2 r 2 



v(t) +uDn(t) 



(B19) 



With D = r/c, and v(t) given by Eq. (|B5 



E vei (r,t) 

Substituting j~ 2 — 1 — /3 2 , and /3 = ujD = ujr/c, 



?E Coul 



euj 2 n(t) 

? 2 ' 

c 7 



(B20) 



E vc i(r, t) — .Ecoui 



-£^Coul = 



euj 2 n(t) 



(B21) 



Using again that -Ecoui = en/R 2 , and that r = v/lj, 



E vcl (r,t) 



'Coul 



= -0 2 



E 



Coul 



= —fiE. 



Coul- 



(B22) 



Therefore representing the electric field at one particle due to the other as merely the instantaneous Coulomb field 
introduces an error only of order /3 4 compared to the exact velocity field. 

The second term on the right hand side of Eq. (IB17[) is the acceleration field term and using Eqs. (|B12I) and (IB13[) 
evaluates to 



(r,t) 



n x ( (n — (3) x j3 



n x ((n — (3) x (3 



e{ojDf 



n x ( (n — /3) x ] 



(B23) 



where here t' = t — 7r/(2w). For the nonrelativistic circular orbit the acceleration is perpendicular to the velocity 
and parallel to n so 
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n x ((n — 0) x 0\ = — n x ( /3 x 



-I3PI3. 



(B24) 



Disregarding initially the delay-related term, it is of interest to determine the relative magnitude of the velocity 
and acceleration field strengths in the system at hand. Consider an electron and positron orbiting in the Bohr model 
ground state with a separation of re = 2ft 2 / (me 2 ) , then 



\E_ 

\E 



acccl | 



vcl 



AS 



j3f3r f3ar v 3 



(B25) 



and 



I $ acccl | 




(?) 




ar 


v 2 


\Bvc\\ 


n 
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c 2 


c 2 



(B26) 



The acceleration fields will be negligible to the analysis in both cases. Corrections to the electric field below order 
f3 2 are not needed. Since the magnetic field enters the dynamics only at the order /3 2 compared to the electric field, 
magnetic field corrections at order j3 2 are also clearly irrelevant here. 

Next, to consider the magnitude of the delay-related terms, and recalling that ljD = j3, the magnitude of the delayed 
term of the acceleration electric field of Eq. (|B23I) can be seen to be a factor of j3 3 smaller than the non-delayed 
component. Therefore the effect of delay on the acceleration fields of the point charges is irrelevant to the present 
analysis. 



The electric acceleration field can thus be approxi- 
mated to order f3 2 as 



-Eaccei(?-,i) 



,,3 1 



e/3 2 



(B27) 



2. Magnetic Field Due to Charge Motion 
From Eq. (|B2[) . the magnetic field "velocity" term is 



B ve i(r,t) 



[n x 0] 



,y2 r 2 L "" " '"Jret ' 



or, using Eqs. (|BT4|) and (|BT6]l . 



(B31) 



B vc i(r,t) 



[(« + uDv) x(/3- uj 2 D 2 n)] . (B32) 



This approximation is accurate to order /3 2 . So, to 
order (3 2 , 



B ve \(r,t) 



n x q 



n-/3 



7 2 (1 - f3 ■ n) 6 R 2 



(B28) 



or 



B vcl (r,t) 



-q(n x (3) 



7 2 (1 - /3 • n) 6 R 2 



(B29) 



As for the case of the electric field described above, for 
the circular orbit (3 ■ n = at the retarded time as well as 
the present and so the magnetic field due to the positron 
charge motion is 



-Bvcl — 



[n X (3 - uj 3 D 3 v x n] 



With ojD = (3 (see above), 

" [(1 + (3 2 )nx(3] 



B 



vel 



B vci = — [(1 - /? 4 )n x 0] 



(B33) 



(B34) 



(B35) 



From Eqs. (|B2| and (|B27|i . the magnetic field "accel- 
eration" term is 



B ve i(r,t) 



-e(n x 0) 



? 2 



(B30) 



B 



acccl 



(r,*) 



n x 



e/3 2 



(B36) 
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Recognizing that the effect of retardation is only an 
order (3 2 change from the field neglecting retardation, 
and that the magnetic acceleration field is overall only 
an order f3 2 contribution to the magnetic field due to 
charge motion, compared to the magnetic velocity field, 
an approximation for the magnetic field that is accurate 
to first order in (3 is 

Bvei = ^ [n x /3] = -jL [v x r] . (B37) 

The magnetic field at the electron due to the orbital 
motion of the positron charge using this approximation 
is then 



B v = -4«„ x r = f i>„ x r p = — ^L p , (B38) 

cr a cr mcr 3 

since r = —2r p = 2r e , and L p = r p x mv p . Noting 
L = 2L„ obtains that 



where the result of Appendix A that the fine structure 
constant a — e 2 /(hc) — 2(3 has been used. The change 
in the electron-positron binding force due to this compo- 
nent of the magnetic field is therefore of order f3 2 , which 
is of the same order as the spin-orbit coupling, and so 
this effect is not negligible to a calculation of the total 
coupling magnitude. However because it is unrelated to 
the angular momentum analysis, due to that it is a ra- 
dial force and so does not contribute a torque, it is of no 
concern to the present analysis. 

3. Magnetic Field Due to Intrinsic Magnetic 
Moments 

The total magnetic field at the electron is due to the 
intrinsic magnetic moment of the positron and to the 
motion of the positron charge. This can be expressed as 
B = B^ + B v . The magnetic field component due to a 
magnetic dipole of moment fi is 



B , 



(B39) 



B,, 



3n (n ■ fi) — /J, 



(B45) 



The torque generated by the magnetic field of the other 
particle charge translational motion will be shown to van- 
ish, but there is a nonvanishing radial force component. 
It is worthwhile to compare this force magnitude to the 
Coulomb force magnitude. The force on the electron 
transiting the magnetic field due to the positron charge 
motion is 



|F| = |-e/3 e x B v \ 
At the Bohr radius 



e 2 f3 



mcr^ 



(B40) 



e 2 /3h 
mere 3 



(B41) 



The relative magnitude of this force to the Coulomb 
attractive force between the electron and positron at the 
Bohr separation is 



-Fcoul 



e 2 /3h 
mere 3 



>'tr 



(B42) 



or 



pn (3h ( 2K 



' Coul 



mere 



(B43) 



or 



LFcoui 2 V he' 



(B44) 



where n = r/r is a unit vector in the direction from 
the dipole to the field point. 

It will be illustrative to consider the relative magnitude 
of the magnetic fields due to orbital motion and intrinsic 
moment, for the quasiclassical Positronium atom at the 
Bohr ground state radius, and for the intrinsic moment 
perpendicular to the orbital plane. Then B^ = —/J-/r 3 
and 



B„ 



\B„ 



\i ges eh 
r 3 2mcr 3 2mcr 3 
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The ratio of the positron intrinsic magnetic field maag- 
nitude to the positron charge velocity magnetic field mag- 
nitude at the electron is 



B„ 



eh 



eh 



B v 2mcr 3 \ mar 3 1 
At the Bohr radius, where L = h, 



h_ 

2L' 



(B47) 



(B48) 



Therefore the intrinsic magnetic field is also not neg- 
ligible to the spin-orbit coupling magnitude. Also, un- 
like the magnetic field due to the other-particle charge 
motion, it can generate a torque for non-perpendicular 
spin orientations relative to the orbital plane, and so has 
already been accounted for in the analysis of angular mo- 
mentum motion. As noted previously, the effect on the 
spin-orbit coupling magnitude of this force in the quasi- 
classical Hydrogen atom was addressed in [l6l - [l9| . 
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4. Electric field due to translational and 
orientational motion of intrinsically-magnetic 
particles 

The vector potential A due to a magnetic moment u 
at a field point outside the source region is [42j 



r 6 



(B49) 



where Tq here is the displacement from the dipole to 
the stationary field point, and r = |r |. Requiring the 
field point to be at the other particle obtains ro = r. 

The electric field is obtained from the vector and scalar 
potentials as 



E 



1 dA 

c~dt 



(B50) 



Although for a stationary magnetic dipole the scalar 
potential vanishes, it cannot be assumed to vanish in the 
present application where the magnetic dipole is trans- 
lating. From Eq. (|B49|) . 



... , dr 1 d 

— =-^(/xxr )- + -3-( / xxr ). 



dA 
~dt 



(B51) 



In the present application only circular orbits are con- 
sidered. For the circular orbit, dr/dt = 0, obtaining that 



dA I f da dr 



(B52) 



Noting that the vector ro here is from the dipole to 
the arbitrary field point, and so v = —dro/dt, the elec- 
tric field component due to the motion of the positron 
magnetic moment relative to the electron, at the elec- 
tron, may be rewritten as 



dA 1 

-fo=^3(t J 'Xr-vxu). (B53) 

The first term will be nonzero for a precessing spin mo- 
ment while the second will be nonzero due to the orbital 
motion of the electron around the proton. The relative 
magnitude of the two terms can be evaluated. The first 
cross product inside the parentheses is upper bounded as 



so, at the Bohr radius 



fir 



< 



7 e s m e 2h 2 
32 h 4 c 2 mc me' 



(B56) 



The second term in the parentheses on the right hand 
side of Eq. (|B52I) can be upper bounded as 



\v x u\ < 



m r 1 / 2 r 1 / 2 mc m 1 / 2 r 1 / 2 mc 
At the Bohr radius 



(B57) 



\v x u\ < 



eh 



-1/2 



to 1 / 2 mc \m r e" 
The ratio of the two term bounds is 



— . (B58) 
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T,Wc~ 2 = k a2 > (B59) 



where a = e 2 /(hc) ~ 1/137 is the fine structure con- 
stant. As derived in the Appendix, v/c for Positronium 
in the Bohr ground state is a/2, so the force due to {x is 
only about of magnitude (v/c) 2 compared to the v x u 
force. Since as shown in Appendix B5 below, the latter 
is already a (v/c) 2 magnitude compared to the Coulomb 
attraction force, the ji related force can be neglected to 
order (v/c) 2 overall. 

To evaluate the electric field term due to the scalar po- 
tential gradient, recall that a translating magnetic dipole 
acquires an electric dipole moment, d, as [35[ 



d = v x u/c, 
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where v here and until stated otherwise below is the 
laboratory frame translational velocity of the magnetic 
dipole (as opposed to the convention used otherwise 
throughout, that v is the electron velocity relative to the 
positron). 

The scalar potential due to an electric dipole is given 

by 
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with (i = es/(mc), with s = \s\. 

An upper bound on s is developed in Appendix D, for 
the quasiclassical Positronium atom model in the Bohr 
ground state, as 
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and 



V$ = — (3(d ■ n)n - d) . 
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where n — r/r is a, unit vector from the source to the 
field point. 

For the electric dipole moment due to the translational 
motion of the magnetic dipole this becomes 
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V<I> = — rr (3((v x fj.) ■ n)n - (v x [i)) . 
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The total electric field due to translational motion 
of the magnetic dipole, accurate to within an error of 
fj 2 times the magnitude of the Coulomb electric field 
strength, is thus 



Evxn = — 3 (3((« x fj,) ■ n)n - 2(v x /*)) . (B64) 



5. Magnitude of Electric Field Due to Motion of 
Intrinsic Magnetic Moments 

It is of interest to determine the relative magnitude 
of the electric field due to translational motion of the 
intrsinsically-magnetic particles to the magnitude of the 
Coulomb electric field. It is needed for assurance that its 
contribution to the total momentum, through the hidden 
momentum that is proportional to the electric field, is 
negligible in the equation of motion of the momentum. 
(It is not negligible through the force that it generates, 
but this is not the only way it enters into the equation of 
motion of the total angular momentum.) 

Eq. (|B64[) provided the total electric field at one 
circular-orbiting particle due to translational motion of 
the other particle magnetic dipole. The ratio this electric 
field magnitude to the Coulomb electric field magnitude 
is bounded as 
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This is of the order of /3 2 . 

Appendix C: Time Rate of Change of the Electric 
Field 

Because the time derivative of the electric field appears 
in the time rate of change of the kinetic momentum (see 
Eq. ([!?)) ), it is desirable to find a compact expression for 
it that is sufficiently accurate. When it was assumed in 
Section IV that it was adequate to simply represent the 
electric field as the Coulomb field in the time-derivative- 
containing terms, it was concluded that these terms con- 
tributed only at order fj 2 and below compaered to the 
translational force term. Therefore, since the needed ac- 
curacy of the time rate of change of the kinetic momen- 
tum is to relative order of (3 2 , it is needed to confirm that 
there are no other contributions to the electric field time 
derivative that are significant to the order of the Coulomb 
field time derivative or higher. This will be confirmed in 
this section. 

Based on the analyses in Appendix B, the electric field 
accurate to order of /3 2 is 
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(The above omits the acceleration field, which is an 
order /3 3 correction, and deviation between the Coulomb 
field and the exact Lienard-Wiechert velocity electric 
field, which is an order /3 4 correction.) 

The time rate of change of the approximate electric 
field is then 
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Assuming a circular orbit obtains that 
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To determine which terms above are relevant to the analysis, consider first that the magnitude of the first term, 
that is the time derivative of the Coulomb field, can be evaluated as 
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This magnitude will serve as a reference for the magnitude bounds of the other terms in (|C3j) 
Next consider 
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The relative magnitude of this term compared to the Coulomb term is thus bounded as 
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At the Bohr radius 
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which is of order ,3 2 here. 
Next consider 
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The relative magnitude of this term compared to the Coulomb term is thus bounded as 
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All other terms in Eq. (|C3j) are of the same magnitude 
as one or the other of these two terms already consid- 
ered. Therefore, approximating the electric field as the 
Coulomb field introduces errors in the time rate of change 
of the electric field at only relative order j3 2 and so is fully 
adequate for the present purpose. 
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To develop an upper bound on s, consider first the 
quantity in the brackets on the right hand side of Eq. 
(|D1[) . Also suppose that the orbital radius is the Bohr 
radius, where L = h. Then 



Appendix D: Upper Bound on the Time Rate of 
Change of the Intrinsic Spin Vector 
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The order of magnitude of the rates of change of the 
electron and positron spin vectors, s e = \s e \ and s p = 
I s p I , is needed at various points in the analysis in order to 
determine the relative significance of various quantities. 
From Eq. ©, 
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Appendix E: Radiation Damping 

Because an electron in a circular orbit will radiate elec- 
tromagnetic energy and momentum, as expected from 
conservation considerations a compensating force known 
as the radiation reaction force will cause the particles to 
decelerate. The radiation reaction is given by [43| 

F ^ = l^- ( E1 ) 

For the circular orbit, for the electron, the radiation 
reaction force magnitude is 




the present analysis, since the angular momentum lost 
to radiation must be reflective of the angular momentum 
change due to the radiation damping forces. Therefore 
the neglect of the radiation fields and radiation damping 
due to charge motion cannot change the central result of 
the present work, that the Thomas precession leads to 
the nonconservation of total angular momentum. 

Radiation damping will also be present associated with 
radiation due to acceleration of the particle magnetic 
dipole moments [44| , and including due to acceleration of 
the electric dipole moment that results when the dipole- 
carrying particles translate. However, these radiation 
damping forces can be shown 45] to be at most of similar 
magnitude to those due to charge motion, at the atomic 
scale. They will also be accompanied by angular mo- 
mentum in the associated radiation fields, equal to the 
amount lost due to the damping forces, and so cannot be 
expected to overcome the nonconservation of total angu- 
lar momentum due to Thomas precession. 
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To compare this with the Coulomb attraction force: 
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So (since a = 2/3 here, by Eq. (|A10|) 'I. the radiation 
reaction force is of order /? 3 smaller than the Coulomb 
attractive force, and so based on other results herein will 
be order /3 smaller than the non-Coulomb forces due to 
the motion of the dipole-carrying particles, and so is neg- 
ligible to the present analysis. This result is consistent 
with the finding that the radiation fields are negligible to 



Appendix F: Effect of Delay on Electric Force 
Induced by Motion of Intrinsically-Magnetic 
Particles 



It has already been shown that representing the exact, 
delayed, electric field at one particle due to the other par- 
ticle charge as simply the instantaneous Coulomb field is 
a very close approximation. It follows directly that the 
effect of delay on the Coulomb attraction is very small. 
Since the instantaneous Coulomb field deviates from the 
retarded velocity electric field only at order /3 4 (see Eq. 
(|B22p ). the ratio of the force error to the overall mag- 
nitude of the Coulomb attraction will also be of order 
j3 , which makes it an order of /3 smaller than the radia- 
tion reaction force as determined relative to the Coulomb 
force by Eq. (|E8|) . Therefore it is of no interest to the 
present analysis. 

A second delay effect can be identified that has not 
yet been addressed, and that is the effect of propaga- 
tion delay on the electric force induced by the motion of 
the magnetic dipoles. Although this latter delay effect 
will also be insignificant to the present analysis, the non- 
radial force it contributes will be seen to be an order of /3 
larger than the non-radial Coulomb force component due 
to delay, at the scale under consideration here. That is, 
it will be about equal to the radiation reaction force, in 
the quasiclassical positronium atom in the Bohr gound 
state. 

For the case of non-relativistic circular orbital motion, 
the propagation delay, D, is approximately that associ- 
ated with the orbital radius as 



D 
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The angular change of the velocity vector during the 
propagation delay is 
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To get an idea of the possible magnitude of the force re- 
sulting from non-instantaneity, we will compute the force 
magnitude based on (|B64[) for the case of a circular orbit 
and spin perpendicular to the orbital plane. The force 
acting on the positron due to the orbital motion of the 
electron intrinsic magnetic moment is then given by 



the separation value at which the two forces will be of 
equal magnitude can be determined. 

Equating the radiation reaction on the electron, as 
given by (|E7|) , with the force due to the moving electron 
intrinsic magnetic moment, that is, evaluating 
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The force magnitude is 
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where the approximation is to neglect the distinction 
between the electron actual and reduced mass. Account- 
ing for propagation delay, the component along the di- 
rection of the proton motion is approximately 
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is the electron-positron separation where the force on 
the positron due to the electron intrinsic magnetic mo- 
ment orbital motion is similar in magnitude to the radi- 
ation reaction force on the electron. For further compar- 
ison, (|F9[) can be expressed in terms of the Bohr radius 



F|| „ sin (1) ^ r -V2 m -m „ l^ r -7/2 m -m 7 

V C' C £C 

(F5) 

or 

„ 1 e\[2 e?\[2h 7/9 , /9 2e 4 h A 9 

F11 « -== — V~ r~ 7/2 m~ 3/2 = — TT-T TTi . (F6) 

c y'mr 2c 2c A 

For comparison, consider the relative magnitude of the 
force due to radiation reaction on the electron. By equat- 
ing the non-radial force component due to delay and the 
radiation reaction force, which will be seen to have dif- 
fering dependencies on the electron-positron separation, 



r a ^r B . (F10) 

So, for the electron spin model and orientation under 
analysis, the separation where the non-radial force on 
the positron is equal to the radiation reaction force on 
the electron is approximately twice the Bohr separation. 
Also, since the non-radial component of the electric force 
increases more rapidly with decreasing range than the 
radiation reaction force, it will be the larger for all smaller 
radii than (|F10|) . 

There is an equal and opposite force on the other par- 
ticle that arises due to the effect of delay on the hidden 
momentum term in Eq. (|32l) . 



[1] E. Rutherford, "The Scattering of a and j3 Particles by 
Matter and the Structure of the Atom," Philososphical 
Magazine, Series 6, Volume 21, May 1911, p. 669-688 

[2] E. Rutherford, "The Structure of the Atom", Philososph- 
ical Magazine, Series 6, Volume 27, March 1914, p. 488- 
498 

[3] Bohr, Nature 92, 231 (1914) 

[4] A. Sommerfeld, "Zur Quantentheorie der Spektrallinien," 

Annalen der Physik 51:1-94 (1916) 
[5] L. de Broglie, "Recherches sur La Theorie des Quanta," 

Ann. de Phys., lOe serie, t. Ill (Janvier-Fevrier 1925). 
[6] G. E. Uhlenbeck and S. Goudsmit, "Spinning Electrons 

and the Structure of Spectra," Nature vol. 117, p. 264- 



265 February 20, 1926 
[7] W. Heisenberg, "ber quantentheoretische Umdeu- 
tung kinematischer und mechanischer Beziehungen," 
Zeitschrift fr Physik, 33, 879-893, 1925 (received July 
29, 1925). 

[8] E. Schrdinger, "An Undulatory Theory of the Mechanics 
of Atoms and Molecules," Physical Review 28 (6) (1926) 
10491070 

[9] Wolfgang Pauli, "Zur Quantenmechanik des magnetis- 
chen Elektrons," Zeitschrift fr Physik (43) (1927) 601-623 
[10] P. A. M. Dirac, "The Quantum Theory of the Electron," 
Proceedings of the Royal Society of London, Series A, 
Containing Papers of a Mathematical and Physical Char- 



34 



acter 117 (778) (1928) 610-624 

[11] L. H. Thomas, "The Kinematics of an Electron with an 
Axis," Philos. Mag. 3, 1-22 (1927) 

[12] D. C. Lush, "Regarding Llewellyn Thomas's paper of 
1927 and the 'hidden momentum' of a magnetic dipole in 
an electric field," |http://arxiv.org/abs/0905.0927| (2009) 
Submitted to Studies in History and Philosophy of Mod- 
ern Physics. 

[13] W. Shockley and R. P. James, " 'Try Simplest Cases' 
Discovery of 'Hidden Momentum' Forces on 'Magnetic 
Currents' ", Phys. Rev. Lett. 18(20), May 1967 (876-879) 

[14] S. Coleman and J. H. van Vleck, "Origin of 'Hidden Mo- 
mentum Forces' on Magnets", Phys. Rev. 171(5), July 
1968 (1370-1375) 

[15] David J. Griffiths, Introduction to Electrodynamics, 
Prentice-Hall, Upper Saddle River, 1999, page 521 

[16] A. L. Kholmetskii, O. V. Missevitch, and T. Yarman, 
"On the classical analysis of spin-orbit coupling in hy- 
drogenlike atoms", Am. J. Phys. 78(4), April 2010 

[17] D. C. Lush, "Comment on 'On the classical analysis 
of spin-orbit coupling in hydrogenlike atoms' ", Am. J. 
Phys. 78(12), December 2010 

[18] A. L. Kholmetskii, O. V. Missevitch, and T. Yarman, 
"Reply to 'Comment on "On the classical analysis of spin- 
orbit coupling in hydrogenlike atoms," by D. C. Lush,' " 
Am. J. Phys. 78(12), December 2010 

[19] D. C. Lush, "On the influence of hidden mo- 
mentum and hidden energy in the classical anal- 
ysis o f spin-orbit coupling in hydrogenlike atoms," 
|http: / /arxiv.org/abs/1009.0495 1 September 2010, Sub- 
mitted to American Journal of Physics 

[20] Lush, 2009, op. cit. See Eq. (77) on page 18. 

[21] Lush, 2009, op. cit. See page 17. 

[22] M. K.-H. Kiessling, "Classical electron theory and con- 
servation laws," Phys. Lett. A 258 (4-6), (1999) pages 
197-204 

[23] A. Schild and J. A. Schlosser, "Fokker Action Principal 
for Particles with Charge, Spin, and Magnetic Moment," 
Jour. Math. Phys. 6(8), Aug 1965, p 1299-1306 

[24] A. Schild and J. A. Schlosser, "Electromagnetic Two- 
Body Problem for Particles with Spin," Jour. Math. 
Phys. 9(6), June 1968, p 913-915 

[25] Thomas, 1927 op. cit. Eq. (4.121). 



[26] V. Bargmann, L. Michel, V. L. Telegdi, " Precession 
of the polarization of particles moving in a homogenous 
electromagnetic field," Phys. Rev. Lett. 2(10), May 1959 

[27] David Hestenes, "The Zitterbewegung Interpretation of 
Quantum Mechanics," Found. Physics. 20(10), (1990) 
1213-1232 

[28] Martin Rivas, Kinematical Theory of Spinning Particles. 
Kluwer Academic, Dordrecht, 2001 (Section 6.1.1, pages 
254-264.) 

[29] J. D. Jackson, Classical Electrodynamics, Third Ed., Wi- 
ley, New York, 1998 

[30] R. H. Good, Jr., "Classical Equations of Motion for a Po- 
larized Particle in an Electromagnetic Field," Phys Rev 
125 (6), (1962) 

[31] V. Hnizdo, "Conservation of linear and angular momen- 
tum and the interaction of a moving charge with a mag- 
netic dipole," Am. J. Phys. 60(3), March 1992 

[32] J. D. Jackson, Classical Electrodynamics, Third Ed., Wi- 
ley, New York, 1998, page 189 

[33] W. Shockley, '"Hidden linear momentum" related to the 
a ■ E term for a Dirac-electron wave packet in an electric 
field,' Phys. Rev. Lett. 20(7), February 1968 (343-346). 

[34] G. Munoz, "Spin-orbit interaction and the Thomas pre- 
cession: A comment on the lab frame point of view," Am. 
J. Phys. 69(5), May 2001 

[35] G. P. Fisher, "The electric dipole moment of a moving 
magnetic dipole," Am. J. Phys. 39, 15281533 

[36] Hnizdo, 1992, op. cit. page 244, Eq. (23) 

[37] G. T. Trammel, "Aharonov-Bohm Paradox," Phys. Rev. 
134(5B), June 1964 

[38] H. Goldstein, Classical Mechanics, Addison- Wesley, 
Reading, 1950, page 133 

[39] Lush, 2009, op. cit. Section IXC, page 20. 

[40] Lush, 2009, op. cit. 

[41] Jackson, 1975, op. cit., page 657 

[42] Jackson, 1998, op. cit., page 186 

[43] Jackson, 1975, op. cit., page 784 

[44] Naoki Itoh, "Radiation reaction due to magnetic dipole 
radiation", Phys. Rev. A 43(2), 15 January 1991, pl002- 
1004 

[45] D. C. Lush, in preparation 



